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now  a  usable  tool  for  many  blast-load  measurement  problems  although  additional  calibration  and  validation 
oqieriments  ve  recommended. 
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SUMMARY 


The  objective  of  this  project  was  to  provide  Ballistic  Research  Laboratory  (BRL)  with  a 
capability  to  measure  explosively  generated  blast  pressure  histories  at  the  surface  of  a  structure 
very  near  the  explosive.  An  example  of  a  situation  of  interest  is  a  1-kg  explosive  charge  spaced 
100  mm  from  a  steel  surface.  The  peak  blast  pressure  on  the  surface  dire^y  under  Uus  charge  is 
about  25  kbar  (2.5  GPa). 

Our  approach  was  to  adapt  the  "flatpack”  stress  gage  to  this  environment.  A  flatpack  stress 
gage  comprises  piezoresistance  foil  elements  encapsulated  in  a  thin  polymer  layer  and  sandwiched 
between  two  steel  strips  welded  at  the  edges.  The  steel  case  allows  the  gage  to  be  welded  to  steel 
armored  vehicles  and  enhances  the  survivability  of  the  gage  elements  under  severe  loading 
conditions.  However,  the  loading  environment  in  this  application  is  so  severe  that  the  gage 
elements  stretch  significantly  and  these  "in-plane  strains"  reduce  gage  accuracy.  In  the  current 
efrbrt,  we  built  flatpack  stress  gages  containing  three  piezoresistance  elements:  a  manganin  (ot 
ytterbium)  element  to  measure  normal  stress  and  two  constantan  elements  to  measure  in-plane 
strains,  liiat  is,  the  nieasured  in-plane  strains  were  used  to  more  accurately  interpret  the  resistance 
history  from  the  manganin  (or  yttrium)  element 

This  approach  required  the  development  of  a  new  flatpack  gage  design  and  the  development 
of  an  analytic^  method  for  interpreting  the  piezoresistance  element  resistance  histories.  The 
flatpack  gages  contain  a  manganin  (or  ytterbium)  element  and  two  constantan  elements  arranged  as 
shown  in  Figure  1.  All  three  sensing  elements  lie  within  a  5-mm  radius.  Figure  2  shows  the 
configuration  of  an  explosive  blast  test  with  a  manganin  flatpack  gage  10.6  cm  from  a  904-g  HE 
charge.  The  resistance  histories  from  this  test  are  plotted  in  Figure  3. 

The  method  for  interpreting  these  resistance  histones  combines  the  piezoresistance  relations 
and  the  constitutive  relations  for  the  three  elements  to  determine  stresses  and  strains  from  the  three 
simultaneous  measurements  of  resistance  change.  The  piezoresistance  and  constitutive  properties 
of  the  element  materials  are  taken  from  previous  work  by  Gupta  and  colleagues.*"^  In  addition  to 
the  requirement  that  the  piezoresistance  properties  and  constimtive  properties  of  each  element  be 
known,  the  attalysis  relies  on  the  assumptions  that  (1)  all  three  elements  experience  the  same 
normal  stress,  (2)  the  elements  are  aligned  in  the  principal  directions,*  and  (3)  the  coupling  of  the 
in-plane  strains  from  the  matrix  (the  insulating  material  inside  the  flatpack)  to  the  elements  is 
known.  The  piezoresistance  relations,  constitutive  relations,  and  coupling  relations  form  a  system 
of  21  scalar  equations.  The  unknowns  are  the  normal  stress  and  in-plane  strains  in  the  matrix  and 
the  six  principal  stresses  and  strains  in  each  element  A  computer  program  named  PIEZOR  was 
written  to  iteratively  solve  these  equations. 

The  data  recorded  in  explosive  blast  experiments  were  analyzed  with  PIEZOR  to  determine 
the  stresses  and  strains  inside  flatpack  gages.  For  the  data  shown  in  Figure  3,  three  analyses  were 
performed  so  that  the  sensitivity  of  the  result  to  two  parameters  could  be  observed.  These 

parameters  are  the  coupling  of  transverse  in-plane  strain  into  the  piezoresistance  elements  [(X33)n] 

and  the  coefficient  relating  resistance  change  to  plastic  strain  (qn)-  solid  lines  in  Figure  4  are 

computed  with  X33  =  0  and  q  =  0  for  all  three  elements.  The  assumption  of  zero  coupling  is 


*This  assumption  is  consistent  with  the  notion  that  the  principal  directions  inside  the  flatpack  are  controlled  by  the 
flatpack  case.  It  does  not  require  that  the  flatpack  itself  be  aligned  with  the  principal  loading  directions. 
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Figure  3.  Resistance  histories  from  a  three-element  flatpack  gage 
in  an  explosive  blast  experiment. 
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Figure  4.  Flatpack  stress  determined  with  PIEZOR. 
Compression  stress  is  negative. 
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consistent  with  the  very  low  (<1%)  coupling  coefficient  for  the  constantan  elements  for  surface- 
mounted  conditions  with  no  normal  stress.^0  Zero  dependence  on  plastic  strain  is  consistent  with 
the  static  properties  of  manganin  and  constantan.  l 

The  dotted  curve  in  Figure  4(a)  is  computed  with  X33  =  1  and  Ti  =  0  for  all  three  elements. 
For  this  case,  the  dependence  of  the  computed  stress  on  the  transverse  strain  coupling  is  small 
because  the  stress  produces  a  large  resistance  change  in  the  manganin  and  the  fla^ack  limits  the 
transverse  strains. 

The  dashed  curve  in  Figure  4(b)  is  computed  with  X33  =  0  for  all  three  elements,  ti  =  0  for 
the  constantan  elements,  but  r\  for  manganin  taken  as  the  function  of  plastic  strain  determined  in 
shock  wave  experiments  (2(X)  ns  duration).^  The  comparison  shows  a  significant  dependence  of 
the  computed  stress  on  T]. 

In  summary,  a  flatpack  stress  gage  containing  piezoresistance  foil  elements  is  v^  well 
suited  for  measuring  blast  pressure  loads  on  the  surface  of  an  armored  vehicle  because  it  can  be 
attached  to  the  target  with  no  penetrations  (by  welding  or  bonding  with  epoxy)  and  it  presents  a 
very-low-profile  obstacle  to  blast  waves.  Two  important  developments  were  used  in  this  project  to 
make  flatpacks  usable  for  this  application.  First,  the  flatpack  design  was  improved  by  using  a 
very-high-strength  steel  case,  incorporating  three  piezoresistance  elements  in  a  single  gage, 
confining  the  piezoresistance  elements  in  a  shallow  cavity  to  mitigate  strains,  and  r^ucing  the  size 
of  the  gage  so  that  all  the  elements  fit  within  a  S-mm  radius.  Second,  an  analytical  procedure  was 
derived  and  a  computer  code  was  written  to  interpret  resistance  histories  from  the  piezoresistance 
elements  in  non-uiuaxial  strain  conditions  based  on  the  constitutive  and  piezoresistance  properties 
of  the  elements.  With  these  developments,  the  flatpack  gage  is  now  a  useful  tool  for  many  blast¬ 
load  measurement  problems. 

A  good  next  step  in  the  advancement  of  the  flatpack  gage  at  BRL  is  to  continue  to  use  it  in 
explosive  blast  environments.  A  repeat  of  the  experiment  shown  here  would  provide  a  measure  of 
reproducibility  of  the  results,  and  a  detailed  calculation  of  this  experiment  would  help  determine  the 
range  of  uncertainty  in  the  mechaiucal  and  piezoresistance  parameters  of  the  elements  (especially 

X33  and  Ti).  Other  experiments  for  which  the  pressures  can  be  calculated  or  measured,  such  as 
contact  explosives  or  explosively  driven  flyer-plate  impacts,  should  aslo  be  performed.  To  further 
the  gage's  utility,  the  flatpack  should  be  tested  under  fragmenting  rounds  and  hardened  as 
necessary  to  perform  in  tiiat  application. 

Additional  basic  work  is  also  needed  to  eliminate  some  of  the  uncertainties  in  the  properties 
of  the  piezoresistance  elements,  and  to  determine  in  more  detail  the  mechanical  response  of  the 
gage  case,  the  insulation,  and  the  elements. 
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SECTION  1 
INTRODUCTION 


The  objective  of  this  project  was  to  provide  BRL  with  a  capability  to  measure  explosively 
generated  blast  pressure  histories  at  the  surface  of  a  structure  very  near  the  explosive.  Au  example 
of  a  situation  of  interest  is  a  1-kg  explosive  charge  spaced  100  mm  from  a  steel  surface.  The  peak 
blast  pressure  on  the  surface  di^tly  under  this  charge  is  about  25  kbar  (2.5  GPa). 

Our  approach  was  to  adt^t  the  "flatpack"  stress  gage  to  this  enviixximent  A  flatpack  stress 
gage  comprises  piezoresistance  foil  elements,  such  as  manganin  or  ytterbium,  and  copper  foil  leads 
encapsulated  in  a  thin  insulating  layer  and  sa^wiched  between  two  steel  strips  weld^  at  the 
edges.  This  package  enhances  the  survivability  of  the  elements  and  leads,  and  it  helps  to  mitigate 
the  in-plane  strains  in  two-  and  three-dimensional  loading  conditions.  Historically,  piezoresistance 
elen^nts  have  been  used  only  in  uniaxial  strain  conditions  (no  in-plane  strain)  for  which  the 
relation  between  normal  stress  and  resistance  change  has  b^n  established  in  gas-gun  impact 
experiments.  However,  the  large  in-plane  strains  produced  in  the  loading  environment  near  an 
explosive  charge  degrade  the  accur^  of  a  flatpack  gage  significantly.  To  overccxne  this  problem, 
we  built  flatpack  stress  gages  containing  constantan  strain  gages  to  measure  the  in-plane  strains, 
and  we  used  this  additional  information  to  more  accurately  interpret  the  resistance  history  from  Ae 
manganin  or  ytterbium  elements. 

This  approach  required  the  development  of  a  new  flaqiack  gage  design  and  the  development 
of  an  analytical  method  for  interpreting  the  piezoresistance  clement  resistance  histories.  In  the  new 
flatpack  gages,  one  ptedommantly  stress-sensitive  element  (manganin  or  ytterbium)  and  two 
predominantly  strain-sensitive  elements  (constantan)  were  arranged  within  a  5-mm  radius.  The 
gage  packages  measured  25  mm  wide,  1.3  mm  thick,  and  about  300  mm  long.  The  new  analytical 
interpretation  method  is  more  than  a  single  conectiem  for  strain.  A  numeric^  algorithm  nam^ 
PIEZOR  was  developed  to  solve  simultaneously  the  combined  set  of  piezoresistance  equations  and 
elastic-plastic  constitutive  equations  for  the  thrM  piezoresistance  elements  in  a  flatpack.  Thus,  in 
addition  to  determining  the  normal  stress,  the  algorithm  provides  all  three  principal  stresses  and 
strains  in  all  three  elements  as  functions  of  time. 

Several  gage-development  experiments  were  perfcHmed  and  two  major  gage-deiix)nstration 
experiments  were  performed.  In  the  demonstration  experiments,  flatpacks  containing  ytterbium 
elements  were  loaded  in  a  gas-gun  impact,  and  flatpaclu  containing  manganin  elements  were 
loaded  with  an  explosive  blast  The  record  from  these  experiments  were  interpreted  with  the 
PIEZOR  algorithin  and  compared  with  calculated  environments.  In  general,  the  three-element 
flatpack  gages  performed  well  and  the  PIEZOR  algmithm  provided  a  straightforward  method  for 
interpreting  the  resistance  histories  for  non-uniaxi^  elastic-plastic  response. 

The  remainder  of  the  repon  is  organized  as  follows:  Sections  2  and  3  describe  the  design 
of  the  flatpack  gages  and  the  dmvation  of  the  analytical  procedure  for  interpreting  their  resistance 
histories.  Sections  4  and  5  describe  results  of  an  explosive  blast  test  and  a  plate  impact  test  cm 
flatpack  gages.  Section  6  presents  our  conclusions  and  offers  recommendations  for  ^ture  work. 
These  sections  are  followed  by  Appendices  A  through  E,  which  describe  the  piezcnesistance 
material  models,  the  PIEZOR  computer  pro^am,  flnite-^ference  computations  of  explosive  blast 
loads,  preliminary  experiments,  and  electronics  omsiderations  fOT  using  flatpack  gages. 
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SECTION  2 
FLATPACK  GAGES 


A  drawing  of  a  three-clement  flatpack  gage  is  shown  in  Figure  5;  photographs  of  flatpacks 
without  covers  are  shown  in  Figure  6.  All  three  sensing  elements  are  roughly  3-mm-square  grids 
and  lie  within  a  5-mm  radius.  The  case  is  made  of  heat-treated  17-7PH  stainless  steel  (Rc  = 
and  comprises  a  base  and  a  cover,  each  measuring  25  mm  wide,  1.3  mm  thick,  and  3(X)  mm  long. 
The  base  has  a  0.36-mm-deep,  20-mm-wide  cavity  into  which  the  elements  and  insulation  are 
placed.  The  elements  and  their  0.05-mm-thick  by  1.3-mm-wide  soldered  copper  leads  are  first 
encapsulated  between  layers  of 0.025-mm-thick  Kapton;  this  sandwich  is  then  potted  in  Shell  Epon 
815  epoxy  in  the  cavity.  Finally,  Ae  base  and  cover  are  welded  together  at  their  edges. 

Specific  assembly  procedures  are  as  follows: 

(1)  Bond  the  elements  to  one  layer  of  Kapton  with  Shell  Epon  815  epoxy. 

(2)  Solder  the  copper  leads  to  the  solder  tabs. 

(3)  Bond  the  other  layer  of  Kapton  over  the  elements  and  leads  with  815. 

(4)  Bond  the  Kapton/element  sandwich  to  the  bottom  of  the  cavity  in  the  17-7PH 
steel  case  with  815,  then  fUl  the  cavity  with  815. 

(5)  the  top  surface  of  the  epoxy  to  make  it  flat  and  flush  with  the  rim  of  the 
steel. 

(6)  Qamp  the  two  steel  pieces  between  copper  chiU  bars.  These  can  be  square- 
section  bars  with  one  hole  down  the  center  and  cold  water  circulating  through 
the  hole.  The  edges  of  the  steel  should  protrude  from  between  the  chill  bars 
by  about  2  mm  for  welding. 

(7)  Weld  the  edges  with  a  TIG  welder.  The  welder  used  in  this  woik  is  a  Miller 
Dialarc  HF-P  with  the  following  settings;  Argon  98%  ,Oxygen  2%,  Straight 
Polarity,  10-55  DC,  Ampere  Adjustment  100,  Renwte  Amperage  Off,  Remote 
Contactor  On,  High  Frequency  Start,  Electrode  Size  15-16.  The  weld  is 
made  with  a  Gullco  Kat  II  tractor  fe^  set  at  10  ipm. 

The  piezoresistance  elements  used  in  the  three-element  flatpack  gages  are  listed  in  Table  1. 
These  elements  are  all  commocially  available,  they  all  have  approxima^y  the  same  grid  geometry, 
and  they  are  all  encapsulated  in  polyimide. 

We  performed  a  three-dimensional  finite-element  analysis  of  this  flatpack  design  in  a  blast 
environnaent,  using  the  DYNA3D  finite-element  code.  *2  In  this  calculation,  the  gage  was  bonded 
to  the  top  surface  of  a  thick  steel  plate,  and  a  surface  pressure  was  applied.  The  spatial  dependence 

of  the  load  was  defined  with  the  Gaussian  distribution  ^cak  =  Poc'^^'^°^  ,  where  Ppcak  is  the  peak 
pressure  at  radius  r,  Pq  =  43  kbar  (4.3  GPa),  and  Tq  =  29mm.  The  pressure  pulse  shape  at  every 
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Table  1.  Piezoresistance  elements  for  flatpack  gages. 


Element 

Model 

Resistance 

(ohms) 

Overall 

Width 

(mm) 

Number 
of  Grid 
Legs 

Grid  Leg 
Length 
(mm) 

Grid  Leg 
Width 
(mm) 

Thickness* 

(mm) 

Manganin 

Micromeasurements 

J2M-SS-110FB-048 

48 

2.8 

12 

3.1 

0.062 

0.005 

Ytterbium 

Dynasen 

Yb4-50-EK 

50 

3.1 

18 

3.1 

0.075 

0.006 

Longitudinal 

Constantan 

Microrrteasurements 

CEA-06-125UN-350 

350 

3.1 

24 

2.8 

0.025 

0.005 

Transverse 

Constantan 

Miaomeasurements 

"CEA-06-125UE-350 

350 

3.1 

24 

2.8 

0.025 

0.005 

The  thickness  is  estimated  from  the  grid  leg  length  and  width,  the  element  resistance,  and  the  nominal 
resistivity  of  the  material. 


radius  was  a  triangle  with  zero  rise  dme  and  a  20-ps  duration.*  The  flatpack  was  oriented  in  the 
radial  direction  and  was  located  with  the  center  of  the  sensing  area  27.8  mm  from  the  center  of  the 
plate.  The  steel  plate,  the  flatpack  case,  and  the  epoxy  were  modeled  as  simple  elastic-plastic 
materials  with  yield  strengths  of  1  GPa,  1.2  GPa,  and  100  MPa,  respectively.  The  piezoresistance 
elements,  leads,  and  Kapton  were  not  modeled  explicitly  because  the  very  small  elements  required 
to  model  them  would  have  gi^tly  increased  confuting  time;  instead,  the  cavity  within  the  flatpack 
was  assumed  to  be  filled  entirely  with  epoxy.  Figure  7  shows  the  normal  stress  and  in-plane 
strains  computed  in  the  epoxy  layer  in  die  flaqiack  gage  at  the  27.8-mm  radius.  At  this  location, 
the  peak  pressure  applied  to  the  surface  of  the  gage  was  17.2  kbar  (1.72  GPa).  The  calculation 

indicates  that  the  gage  responds  in  a  damped  oscillatory  mode  with  a  period  of  about  5  ps.  This 
oscillation  was  also  observed  in  the  blast  experiments,  described  in  Section  4  and  Appendix  D. 

The  calculated  in-planc  strains  are  less  than  1%,  also  consistent  with  experimental  observation. 

The  oscillation  in  the  gage  stress  is  caused  by  the  motion  of  the  cover  plate  on  the  epoxy- 
filled  cavity.  The  formula  for  the  period  of  this  motion  is 


where  p  and  h  are  the  density  and  thickness  of  the  cover  plate,  and  K,  G,  and  d  are  the  bulk 
modulus,  shear  modulus,  and  depth  of  the  epoxy  layer.  (Using  p  =  7.9  g/cm^,  h  =  1.27  mm,  K  = 

3.9  GPa,5  G  =  1.5  GPa,^  and  d  =  0.36  mm,  we  calculate  the  period  as  4.9  ps).  The  oscillation  of 
the  gage  stress  is  accentuated  in  the  finite-element  calculation  for  two  reasons.  First,  the  rise  time 

*This  pressure  loading  was  based  on  the  results  of  a  nrule-difTerence  computation  of  a  16S-g  C-4  explosive  charge 
over  a  steel  plate  (see  Appendix  Q. 
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Figure  7.  Three-dimensional  finite-element  analysis  of  a  flatpack  gage 
under  blast  loading. 
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of  the  blast  pressure  in  the  calculation  was  zero,  whereas  an  actual  blast  pressure  loading  has  a 
flnite  rise  time.  Second,  this  gage  is  designed  for  use  in  larger  scale  blast  environments  for  which 
the  duration  will  be  a  factor  of  2  or  more  longer  than  in  the  calculation,  which  corresponds  to  our 
sub-scale  developmental  experiments. 

In  theory,  the  magnitude  of  stress  oscillations  could  be  reduced  by  increasing  the  natural 
frequency  of  the  gage.  Cfee  way  to  accomplish  this  would  be  to  decrease  the  cover  thickness  and 
cavity  depth  (h  and  d)  and  to  increase  the  stiffness  of  the  insulating  layers  (K  and  G).  Other 

options  would  be  to  decrease  the  density  of  the  cover  plate  (p),  by  replacing  the  high-strength 
stainless  steel  cover  with  titanium,  or  to  miniaturize  the  entire  package. 

In  an  attempt  to  reduce  the  magnitude  of  the  stress  record  oscillations,  we  built  two  gages 
with  the  cavity  depth  cut  in  half  to  0.18  mm  (0.007  inch).  In  addition,  the  Kapton  insulation  was 
replaced  with  mica,  which  has  an  elastic  modulus  similar  to  that  of  steel.  In  one  gage,  the  cover 
thickness  was  also  reduced  by  a  factor  of  2  to  0.64  mm  (0.025  inch).  In  an  explosive  blast 
experiment  (Test  6,  Section  4)  the  mica-insulated  gage  with  the  standard  cover  showed  less 
overshoot  than  a  Kapton-insulated  gage,  but  the  in-plane  strains  were  larger  and  the  gage  failed 
before  the  measurement  was  complete.  The  larger  strains  were  not  expected,  and  this  may  not  be  a 
general  result.  However,  mica  is  a  layered,  anisotropic  mineral  and  it  may  allow  more  shear 
deformation  between  layers  than  assumed  in  this  design.  The  gage  with  the  thin  cover  failed  at  the 
time  of  arrival  of  the  blast 
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SECTION  3 


PIEZORESISTANCE  ANALYSIS 


Because  some  in-piane  strain  in  the  manganin  and  3merbium  elements  is  inevitable,  a 
method  is  needed  to  account  for  these  strains  when  interpreting  the  resistance  histories  to  determine 
stresses.  Furthermore,  because  the  elements  will  undergo  plastic  deformation,  the  analysis  of 
piezoresistance  gage  records  must  include  elastic-plastic  element  response.  This  section  describes 
the  PIETOR  algorithm,  an  analytical  procedure  developed  in  this  project  to  interpret  piezcxesistance 
flatpack  gage  records. 

As  described  earlier,  the  measurement  strategy  adopted  here  is  to  record  resistance  histories 
from  three  independent  piezoresistance  elements  positioned  at  the  same  "point"  within  the  flatpack 
and  oriented  in  the  principal  directions  of  the  insulation  matrix  stress  and  strain  fields.  In  the 
example  shown  in  Figure  8,  one  manganin  element  and  two  constantan  elements  oriented  at  right 
angles  to  one  another  are  used  to  measure  normal  stress  and  two  in-plane  strains.  The  manganin 
element  is  primarily  sensitive  to  stress  normal  to  the  element;  the  constantan  elements  are  primarily 
sensitive  to  strain  in  the  direction  of  the  element  grid.  Each  element,  however,  has  sensitivity  to 
three  principal  stresses  and  strains.  The  analysis  simultaneously  solves  the  piezoresistance 
equations  and  the  elastic-plastic  constitutive  equations  for  the  three  elements  to  determine  stresses 
and  strains  from  three  simultaneous  measurements  of  resistance  change. 

In  addition  to  the  requirement  that  the  piezoresistance  properties  and  constimtive  properties 
of  each  element  be  known,  ^e  analysis  relies  on  the  assumptions  that  (1)  all  three  elements 
experience  the  same  normal  stress,  (2)  the  elements  are  aligned  in  the  principal  directions,  and 
(3)  the  coupling  of  the  in-plane  strains  from  the  matrix  (the  insulating  material  inside  the  flatpack) 
to  the  elements  is  known.  The  piezoresistance  and  constitutive  properties  of  the  element  materials 
are  taken  from  previous  work  by  Gupta  and  colleagues.  The  first  assumption  is  supported  by 
the  results  of  finite-element  calculations  performed  by  Ito  and  Muki.^  The  second  assumption  is 
consistent  with  the  notion  that  the  principal  directions  inside  the  flatpack  are  normal  to  the  planes  of 
symmetty  of  the  case,  i.e.,  longitudinal,  transverse,  and  through-the-thickness.  The  third 
assumption  is  a  statement  of  the  so-called  inclusion  problem,  for  which  a  general  solution  is  not 
available.  For  strains  in  the  direction  of  the  element  grid  legs,  one-to-one  coupling  is  safely 
assumed  because  of  the  tremendous  aspect  ratio  of  the  ^d  legs.  Although  the  coupling  of 
transverse  strain  to  the  element  grids  is  much  less  certain,  the  analysis  suggests  that  the  resolved 
stress  is  not  terribly  sensitive  to  this  relation. 

As  shown  in  Figure  8,  an  x,  y,  z  coordinate  system  is  fixed  in  the  flatpack,  and  a  1, 2, 3 
coordinate  system  fixed  in  each  element  with  the  1  direction  parallel  to  the  element  grid,  the  2 
direction  transverse  to  the  element  grid,  and  the  3  direction  normal  to  the  element  grid.  The 
increments  of  principal  stress  and  strain  in  the  nth  element  are  grouped  into  two  one-dimensional 
arrays  Sn  and  tn>  and  a  one-dimensional  array  f  contains  the  increments  in  normal  stress  and  in¬ 
plane  strains  in  the  matrix. 
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Figure  8.  Orientation  of  three  piezoresistance  elements  in  a  plane. 


The  rotation  matrix  On  relates  the  matrix  stresses  and  strains  in  the  flatpack  coordinate 
system  to  those  in  the  n-th  element  coordinate  system.  From  Mohr’s  circle  for  two-dimensicmal 
strain,!  1 


Qn  = 


■  1 
0 
0 


0 

cos^Gn 

sin^Gn 


0 

sin^Gn 

cos^Gn 


(2) 


where  Gn  is  the  angle  between  the  X-axis  and  the  In-axis. 

The  coupling  matrix  An  relates  the  normal  stress  and  in-plane  strain  increments  in  the  nth 
element  to  those  in  the  matrix. 


An  = 


1  0  0 
0  1  0 

Lo  0  X33  J 


(3) 


Thus,  the  normal  stress  and  in-plane  strain  increments  in  the  nth  element  are  related  m  those  in  the 
matrix  by 


Sn  =  An  On  f 


(4) 


'Rie  form  of  An  implies  that  the  normal  stress  and  in-plane  strain  increments  along  the 
length  of  the  element  grid  are  identical  to  the  matrix  stress  and  strain  increments  in  those  directions, 
but  the  transverse  in-plane  strain  increment  in  the  element  is  coupled  to  the  matrix  strain  increment 

by  the  faaor  X33.  The  An  matrix  could  be  made  more  general  to  represent  an  inclusion  analysis  of 
the  element  in  the  matrix.  (By  setting  X22  =  ^33  =  0,  we  impose  uniaxial  strain  conditions.) 

The  incremental  elastic-plastic  constitutive  relation  for  the  nth  element  is  (with  the  subscript 
n  omitted  for  clarity 

Aai  =(K-»-|gN)  Aei-KK-|gN)  Ae2  +  (K-|GN)  Ae3  +  aIo(N- 1) 

Ao2  =  (K  -  |gN)  Aei  (K  +  |gN)  Ae2  +  (K  -  |gN)  A£3  +  aio(N  - 1) 

Aa3  =  (K  -  |gN)  Aei  +  (K  -  |gN)  Ae2  +  (K  |gN)  Ae3  +  a5o(N  -  1)  (5) 

where  Aoi  and  Aei  are  principal  stress  and  strain  increments,  aio  and  0(0  are  initial  values  of  the 
principal  stresses  and  stress  deviators,  and  K  and  G  are  the  bulk  and  shear  moduli.  N  is  given  by 
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(H  +  3G  +  3G),  where  H  is  the  hardening  modulus,  5©  is  the  initial  value  of  the  effective 

oN 

stress  expressed  as  +  (olo)^  +  (ojo)^l .  and  is  the  effective  stress  computed 

with  entirely  elastic  strain  increments,  for  which  the  stress  deviators  are  =  aio  +  2GA£i. 

Using  matrix  notation  and  rearranging  terms,  we  can  rewrite  the  constitutive  relation 


Aoi 

'oio' 

Aa2 

-(N-1) 

aj>o 

_ACT3_ 

-CT^O- 

mil  mi2  mi3- 

m2i  TD22  r.  ^ 

msi  m32  m33 . 


Aei 

Ae2 

_Ae3_ 


(6) 


and  make  the  independent  quantities  increments  in  normal  stress  and  in-plane  strains 


Aoi 

■(N  -  l)alo' 

(N  -  l)aio 

Ao2 

=  M  1  Aei 

(N  -  l)aio 

-M* 

0 

_Ae3_ 

LAe2. 

L  0 

0 

where 


r  (™) 

^33^ 

“  m33  ' 

^  m33  ^ 

M*  = 

(^) 

^m33^ 

(m2i-“^) 

ni33 

ni22m22) 

m33  ^ 

(-^) 

(.mil) 

"133 

(- 

m33"  -J 

For  the  nth  element,  this  equation  condenses  to 


tn  =  Mn  Sn+  bn 


where 


‘(N-  Dole' 

(N  -  1)0^0 

(N  -  l)aio 

-  Mn 

0 

L  0 

n 

0 

(7) 


(8) 


(9) 


(10) 
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The  piezoiesistance  equation  that  relates  the  incremental  change  in  resistance  of  the  nth 

element  ARn,  nonnalized  by  the  initial  value  of  resistance  (Ro)n>  to  an  increment  in  principal 
stresses  and  strains  is 


ARn 

(Ro)n 


Vs  D 

=  ai  Aasn  +  Aein  -  Ae2n  +  (“n  +  2pn)Aoin  +  an  Aa2n  “  A£3n  +  ^In  -2"  (1 1) 


where  an,  Pn,  and  "Hn  are  piezoresistance  properties  of  the  nth  element  and  Ae^  is  the  increment  in 
the  plastic  strain  in  the  nth  element  This  relation  can  be  written  in  the  form 


"Aasn" 

+  [  (an  +  2Pn)  an  -  1  ] 

'Aain 

rn  -  [  ttn  1  -1  ] 

Aein 

Aa2n 

.Ae2n- 

-Aesn- 

+  ’i"T 


(12) 


or 


fn  =  Pn  Sn  +  Qn  tn +  ‘nn^  (13) 


Equations  (4),  (9),  and  (13)  form  a  system  of  21  scalar  equations.  The  unknowns  are  the 
normal  stress  and  in-plane  strains  in  the  matrix,  contained  in  f,  and  the  6  principal  stresses  and 
strains  in  each  element,  contained  in  Sn  and  tn-  These  are  solved  by  substituting  (4)  and  (9)  into 
(13)  to  get 


rn  =  Pn  An  fin  f  +  Qn  Mn  An  fin  f  +  Qnbn+  Tin  “2“  ^^n 
which  can  be  rewritten  as 

rn  ~  “Hn  ^  Ae^  -  Qn^n  =  fPn  An  fin  +  Qn  Mn  An  fin!  ^  (15) 


or,  condensing  further 


r*  =  Wf 


(16) 


yj  n 

where  r*  is  a  3  X  1  array  whose  nth  element  is  rn  -  ‘Hn'T'  ~  Qn^n  and  W  is  a  3  x  3  array 

^  n 

whose  nth  row  is  given  by  Pq  Aq  fin  +  Qn  Mn  An  fin-  Equation  (16)  is  solved  iteratively  for  f; 
then  equations  (4)  and  (9)  provide  Sn  and  tn.  A  computer  program  named  PIEZOR  was  written  to 
perform  this  task. 
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The  PBEZOR  progr^  was  tested  with  a  sample  problem  generated  with  a  NIKE2D  finite- 
element  computation.  In  this  problem,  constantan  and  manganin  elements  aligned  like  those  in 
Figure  8  were  subjected  to  the  bounda^  conditions  plotted  in  Figures  9(a)  and  (b).  The  normal 
stress  ramped  linearly  to  10  kbar  (1  GPa)  and  remained  at  that  level;  the  longitudinal  strain 

oscillated  in  a  saw-tooth  wave  form  with  an  amplitude  of  ±1%;  the  transverse  strain  remained  zero. 
Because  these  boundary  conditions  were  applied  directly  to  each  element  (not  through  a  matrix 

material),  An  is  the  identity  matrix  for  each  element.  The  stresses  and  strains  computed  with 
NIKE2D  were  used  in  equation  (13)  to  obtain  resistance  histories  for  each  element  These  are 
shown  in  Fi^e  9(c).  Note  that  for  these  conditions  the  output  of  the  manganin  gage  is  far  from 
constant  during  the  phase  of  constant  stress. 

The  resistance  histories  were  input  to  PIE2X)R,  which  calculated  the  element  stresses  and 
strains.  Comparisons  of  the  output  from  PIEZOR  with  the  output  from  NIKE2D  are  shown  in 
Figures  10  through  13.  All  the  stresses  and  strains  in  each  element  were  computed  with  very  good 
accuracy  (<1%  error). 
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Rgure  9.  Applied  boundary  conditions  and  computed  resistance 
histories  for  sample  piezoresistance  analysis  problem 
(T|  -  0). 
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Figure  10.  Manganin  element  (Ml)  stresses  computed  with  NIKE2D 
finite  element  code  (solid)  and  PIEZOR  piezoresistance 
analysis  (dashed). 
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Figure  1 1 .  Manganin  element  (Ml)  strains  computed  with  NIKE2D 
finite-element  code  (solid)  and  PIEZOR  piezoresistance 
analysis  (dashed). 
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Figure  12.  Constantan  element  (Cl)  stresses  computed  with  NIKE2D 
finite  element  code  (solid)  and  PIEZOR  piezoresistance 
analysis  (dashed). 
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Figure  1 3.  Constantan  element  (Cl )  strains  computed  with  NIKE20 
finite  element  code  (soKd}  arxi  PIEZOR  piezoresistance 
analysis  (dashed). 


SECTION  4 

CLOSE-IN  BLAST  EXPERIMENT 


Four  flatpack  gages  of  three  different  designs  were  used  in  the  close-in  blast  experiment 
(Test  6).  Each  flatpack  contained  a  manganin  element  and  two  constantan  elements  in  an  epoxy- 
fllled  cavity.  The  four  flatpack  gages  were  welded  to  the  top  surface  of  a  25-cm-diameter  steel 
cylinder  in  the  positions  shown  in  Figure  14.  Flatpacks  FP-1,  FP-2,  and  FP-3  were  positioned  at 
a  2.S4-cm  radius  (measured  to  the  center  of  the  manganin  element);  flatpack  FP-4  was  positioned 
at  a  S.08-cm  radius. 

The  differences  in  the  designs  were  in  the  insulating  material,  the  depth  of  the  cavity,  and 
the  thickness  of  the  stainless  steel  cover  plate.  In  two  of  the  flatpacks  (FP-1  and  FP-4)  the 
elements  were  insulated  with  Kapton  film  and  placed  in  a  0.35-mm-deep  cavity  (0.014  in.)  filled 
with  epoxy.  These  flatpacks  both  used  1.27-mm-thick  cover  plates.  In  the  other  two  flatpacks 
(FP-2  and  FP-3),  the  elements  were  insulated  with  mica  film  and  placed  in  a  0.18-mm-deep  cavity 
(0.007  in)  filled  with  epoxy.  One  of  these  (FP-2)  also  used  a  1.27-mm-thick  cover  plate;  Ae  other 
(FP-3)  used  0.64-mm-thick  cover  plate. 

Figure  15  shows  the  configuration  of  the  explosive  blast  test  The  904-gram  explosive 
charge  was  an  8.8-cm-diameter  by  8.8-cm-tall  cylinder  of  Composition  B,  provided  by  BRL  (ID 
No.  ABY  90A0E1S027).  The  charge  was  positioned  over  the  25-cm-diameter  steel  cylinder  at  a 
standoff  of  15  cm,  measured  to  the  center  of  gravity  of  the  charge  (10.6-cm  to  the  bottran  of  the 
charge).  The  charge  was  supponed  by  a  2.54-cm-tiiick  disk  of  styrofoam  having  a  6-cm-diameter 
centi^  hole  and  an  8.8-cm-diameter,  1.27-cm-deep  counter  bore.  Detonation  was  initiated  with  a 
Reynolds  RP-1  detonator  placed  at  the  top  center  point  of  the  charge.  This  configuration  was 
chosen  to  match  a  test  perfonned  on  an  instrumented  target  plate  at  FMC  Corporation.^  The 
nominal  peak  pressure  on  the  target  directly  undo*  the  charge  was  expected  to  be  about  27  kbar 
(2.7  GPa),  based  on  calculations  performed  at  Southwest  Research  Institute  (SWRI).^ 

The  histories  of  resistance  change  were  recorded  on  digital  oscilloscopes  at  a  sampling  rate 
of  10  points  per  ps,  which  provided  a  recording  window  extending  to  about  190  ps  after  time  of 
arrival  (TOA).  FP-1  provided  signals  for  about  100  ps  after  TOA,  FP-2  provided  signals  for 

about  20  ps  after  TOA,  FP-3  failed  within  5  ps  after  TOA,  and  FP-4  provided  signals  for  the  full 
duration  of  the  recording  window.  The  resistance  histories  from  the  elements  in  FP-1,  FP-2,  and 
FP-4  are  plotted  in  Figures  16  through  18.  In  the  time  scale  of  these  plots,  the  detonator  was 
triggered  at  t  =  0. 

The  data  from  FP-1,  FP-2,  and  FP-4  were  analyzed  with  the  PIEZOR  computer  program  to 
determine  the  actual  stresses  and  strains  inside  the  flatpack  gages.  Recall  that  in  the  PIEZOR 
analysis  an  assumption  must  be  made  regarding  the  coupling  of  transverse  in-plane  strain  into  the 

piezoresistance  elements  (see  Section  3  for  the  definition  of  the  A  matrix).  A  very  low  (<1%) 
coupling  coefficient  is  quoted  by  Micromeasurements  for  the  strain  gage  elements,  but  this 
coefficient  is  determined  only  for  surface-mounted  gages  with  no  non^  stress.  The  solid  curves 

plotted  in  Figures  19  through  23  use  the  assumption  that  transverse  strain  is  fully  coupled  (X33  = 

1);  the  dashed  curves  use  the  assumption  that  transverse  strain  is  completely  uncoupled  (X33  =  0). 
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Figure  1 4.  Fla4>ack  position  in  Test  6. 
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Figure  1 5.  Configuration  of  Test  6. 


TIME  (4S) 

(c)  Constantan  2 


Figure  16.  ResistarKe  histories  from  FP-i 
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Figure  18.  Resistance  histories  from  FP-4. 
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Figure  19.  Normal  stresses  in  FP-1  and  FP-2. 
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Figure  21 .  In-piane  strains  in  FP-1 . 
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(a)  LongitucSnal  strain 


These  two  assumptions  are  approximate  bounds  on  the  actual  coupling.  In  these  analyses  the 
piezoresistance  hardening  coefficient  t\  for  manganin  was  assumed  to  be  a  function  of  the  plastic 
strain,  as  determined  in  shock-loading  experiments  (200  ns  duration).'*  See  Appendix  A  for  the 
material  properties  and  piezoresistance  coefficients. 

Several  observations  can  be  made  from  the  data.  The  resolved  late-time  stresses  in  FP-1 
[Figure  19(a)]  and  the  early-time  stresses  in  FP-4  [Figure  20]  appear  to  suppon  the  assumption  of 

zero  coupling  of  transverse  strain.  That  is,  the  dashed  curves  in  the  figures  (X33  =  0)  are  more 

physically  plausible  because  we  do  not  expect  the  stress  at  2.54  cm  to  become  negative  after  30  ps 
and  we  do  expect  a  peak  stress  of  about  4  kbar  (0.4  GPa)  at  5.08  cm.  For  both  FP-1  and  FP-2,  in 
which  the  stress  produces  a  large  resistance  change  in  die  manganin,  the  dependence  of  the 
computed  stress  on  the  transverse  strain  coupling  assun^tion  is  small  because  the  flatpack  limits 
the  transverse  strains  sufficiently.  In  the  case  of  FP-4  this  is  not  true,  and  it  suggests  that  a 
ytterbium  element  (more  sensitive  to  stress)  would  have  been  a  better  choice  for  the  5.08-cm  range. 

Figure  19  shows  that  the  stresses  in  FP-1  and  FP-2  (up  to  the  time  of  failure  of  FP-2)  are 
very  nearly  the  same,  which  is  consistent  with  the  symmetry  of  the  experiment.  The  peak  stress  in 
FP-4  (Figure  20)  is  much  lower  than  those  in  FP-1  and  FP-2,  which  is  consistent  with  the  rapid 
decay  of  the  blast  wave  with  radius.  Figures  21  and  22  show  that  during  the  time  of  interest  both 
the  longitudinal  and  transverse  strains  in  FP-1  are  about  half  those  in  FP-2,  suggesting  that  the  FP- 
1  design  is  the  better  of  the  two.  (FP-1  also  survived  much  longer.) 

One  more  analysis  was  performed  with  the  data  from  FP-1.  In  this  analysis,  the 
piezoresistance  hardening  coefficient  r\  for  manganin  was  assumed  to  be  zero,  which  is  consistent 
with  data  from  quasistatic  tests.  ^  Figure  24  shows  that  this  assumption  raises  the  resolved  stress 
in  the  flatpack  at  late  time  to  about  2  kbar  (0.2  GI^).  This  result  suggests  that  T)  lies  between  the 
values  detemtined  from  static  experiments  ^  and  gas-gun  experiments.  3 
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Figure  24.  Normal  stress  in  FP-1 . 
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SECTION  5 

PLATE-IMPACT  EXPERIMENT 

The  objective  of  the  plate-impact  test  was  to  demonstrate  the  performance  of  the  flatpack 
gage  and  the  PIEZOR  algorithm  under  conditions  for  which  the  applied  loading  is  known. 
Symmetric  plate-impact  experiments  provide  a  known  value  of  incident  stress  if  the  flyer  and  target 

density,  the  flyer  velocity,  and  the  shock  velocity  are  measured  (o  =  ^  pcv).  Uniaxial  strain 

loading  is  ntaintained  until  relief  waves  from  the  perimeter  of  the  flyer  and  target  arrive  at  a 
measurement  location.  Although  an  ideal  loading  for  a  blast-pressure  gage  would  be  a  known 
high-pressure  shock  wave  in  a  gas,  plate-impact  experiments  use  solid  or  liqtud  media.  Solid 
flyers  and  targets  are  by  far  easier  to  use.  Thus,  our  approach  was  to  place  gages  on  the  front  of  a 
steel  target,  cover  them  with  a  PMMA  buffer  plate,  and  impact  this  assembly  with  a  PMMA  flyer. 

The  design  of  the  test  is  shown  in  Figure  25.  The  25.4-mm-thick  flyer  plate  produces  a 
pulse  duration  of  about  20  ps.  The  25.4-mm-thick  buffer  plate  delays  the  arrival  of  the  rarefaction 

from  the  upstream  face  of  the  flyer  at  the  gage  plane  by  another  20  ps.  The  thickness  of  the  steel 
target  plate  is  such  that  the  rarefaction  from  the  rear  of  the  target  arrives  at  the  gage  plane  at  about 
the  same  time.  The  300-mm-diameter  of  the  flyer  and  target  was  selected  to  allow  uniaxial  strain 

conditions  to  maintain  for  about  30  ps  Gonger  than  the  incident  pulse  duration).  Two  flatpack 
gages  (FP-1  and  FP-4)  and  two  Kapton-covered  element  clusters  (FP-2  and  FP-3)  were  mounted 
on  the  steel  target  plate.  Shallow  channels  were  machined  into  the  surface  of  the  PMMA  buffer 
plate  to  accommodate  the  gages.  These  were  filled  with  Epon  815  epoxy  when  the  plates  were 

bonded  together.  The  measured  impact  velocity  was  44.7  m/s.  The  measured  tilt  was  9  ps  from 
top  to  bottom. 

One-,  two-,  and  three-dimensional  pretest  computations  were  performed  (see  Appendix  F). 
These  predicted  that  the  stress  inside  the  flatpacks  woidd  be  only  slightly  perturb^  from  the 
reflected  stress  pulse  at  the  steel  target  face.  More  specifically,  the  registration  of  stress  into  the 
flatpack  was  predicted  to  be  one-to-one,  and  the  oscillations  in  the  stress  caused  by  the  motion  of 

the  cover  plate  were  predicted  to  be  about  ±0.2  kbar  {20  MPa).  For  the  measured  impact  velocity, 
the  predicted  reflect^  stress  was  1.35  kbar  (135  MPa). 

As  shown  in  Figure  26,*  the  resistance  histories  from  the  ytterbium  elements  inside  the 
flatpacks  were  very  different  fiiom  those  mounted  directly  on  the  steel  plate,  and  neither  set 
matched  the  predictions.  First,  the  bare  (Kapton  covered)  element  records  indicate  a  very  much 

more  rounded  wave  front  than  expected,  although  the  pulse  duration  (about  17.4  ps)  and  the 
arrival  time  of  the  second  rarefaction  are  about  as  predicted.  The  elements  inside  the  flatpacks 
indicate  a  much  lower  equilibrium  level  and  a  large  overshoot .  In  Figure  27,  these  resistance 
histories  are  converted  to  stress  using  the  assumption  of  uniaxial  strain  and  element  properties 
based  on  Gupta's  results  for  foils**  (see  Appendix  A).  This  interpretation  indicates  that  the 
equilibrium  stress  level  inside  the  flatpacks  was  about  half  of  that  at  the  steel  target  face. 


*Zero  time  in  the  plots  corresponds  to  first  sense  of  impact  (9-ps  tilt). 

**The  ytterbium  elements  are  vapor-deposited  rather  than  foils,  so  the  properties  may  be  different 
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Rgure  27.  Stresses  computed  for  uniaxial  strain  conditions. 
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The  measured  tilt  was  9  (top  to  bottom).  Tilt  may  be  the  reason  for  non-simulaneity  of 
the  signals  from  the  four  gages,  but  tilt  does  not  account  for  the  nature  of  the  records. 

The  constantan  element  records  are  shown  in  Figures  28  through  31.  For  the  bare  gages, 
these  signals  look  suspiciously  like  electromagnetically  induced  signals  because  during  the  uniaxial 

strain  phase  (before  about  30  ps)  the  constantan  output  resembles  the  time-derivative  of  the 

ytterbium  output  (except  for  the  2.5-|is-late  time-of-arrival  for  the  FP-3  ytterbium  element,  which 
cannot  be  explained).  The  large  excursions  in  the  records  from  the  bare  constantan  elements  at 

about  45  [IS  coincide  with  the  arrival  of  the  rarefaction  from  the  rear  of  the  target.  The  constantan 
elements  inside  the  flatpacks  indicated  relatively  small  strains  with  some  oscillations  during  the 
uniaxial  strain  phase.  In  contrast  to  the  bare  elements,  the  flatpack  case  mitigated  the  late-time 
strains. 


The  data  from  flatpacks  FP-1  and  FP-4  were  interpreted  with  PIEZOR  three  ways.  First, 
the  transverse  strain  coupling  (X33)  was  assumed  to  be  zero  and  the  piezoresistance  strain¬ 
hardening  coefficient  (il)  was  set  to  zero.  Second,  X33  =  1  with  Ti  =  0.  Third,  >.33  =  0  with  ti  a 
function  of  the  plastic  strain  (see  Appendix  A).  These  interpretations  are  shown  in  Figures  32  and 
33.  None  of  these  interpretations  brings  the  flatpack  stresses  into  agreement  with  the  free-field. 
The  data  from  FP-2  and  FP-3  were  not  analyzed  with  PIEZOR  because  the  measured  strains  are 
zero  during  the  uniaxial  strain  phase  (for  which  the  stress  is  known). 

We  do  not  fully  understand  these  results.  It  may  be  that  the  compliance  of  the  flarpack  is 
much  greater  than  assumed  in  the  pretest  calculation.  That  is,  because  of  this  compliance,  the  full 
stress  in  the  PMMA  buffer  was  not  transmitted  to  the  sensing  region  of  the  flatpack.  Also,  in  our 
calculation  we  assumed  that  the  PMMA  and  epoxy  remained  elastic  for  this  loa^ng  (sec  Appendix 
F  for  details).  Plastic  response,  especially  inside  the  flatpack,  may  be  a  source  of  lower 
compliance.  A  set  of  posttest  calculations  would  likely  provide  a  much  better  understanding  of  this 
test.  However,  these  calculations  could  not  be  includ^  in  the  current  contract. 
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Figure  28.  Constantan  element  resistance  histories  for  FP-2  (Kapton-covered). 
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Figure  30.  Constantan  element  resistance  histories  for  FP-1  (flatpack  gage). 
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Figure  31 .  Constantan  element  resistance  histories  for  FP-4  (flatpack  gage). 
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Figure  32.  PIEZOR  interpretation  of  FP-1  records. 
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(a)  Effect  of  transverse  strain  coupling  with  -i\ «  0. 


(b)  Effect  of  plastic  strain  with  ^33  «  0. 
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Figure  33.  PIEZOR  interpretation  of  FP-4  records. 
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SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


A  flatpack  stress  gage  containing  piezoresistance  foil  elements  is  very  well  suited  for 
measuring  blast  pressure  loads  on  the  siMace  of  an  armored  vehicle  because  it  can  be  attached  to 
the  target  with  no  penetrations  (by  welding  or  bonding  with  epoxy)  and  it  presents  a  very-low- 
profile  obstacle  to  blast  waves.  However,  at  the  beginning  of  this  project,  the  usefulness  of 
flatpack  gages  for  this  application  was  limited  primly  b^use  large  in-plane  strains  in  the 
elements  would  significantly  perturb  the  element  output  and  prevent  it  from  being  properly 
interpreted. 

In  this  project  we  accomplished  two  major  developments  to  make  flatpacks  usable  for  this 
application.  First,  the  flatpack  design  was  improved  by  using  a  very-high-strength  steel  case, 
incorporating  three  piezoresistance  elements  in  a  single  gage,  confining  die  piezoresistance 
elements  in  a  shallow  cavity  to  mitigate  strains,  and  reducing  the  size  of  the  gage  so  that  all  the 
elements  fit  within  a  5-mm  radius.  Second,  an  analytical  p^edure  was  deriv^  and  a  ctxnputer 
code  was  written  to  interpret  resistance  histories  from  the  piezoresistance  elements  under 
nonuniaxial  strain  conditions  based  on  the  constitutive  and  piezoresistance  properties  of  the 
elements.  With  these  developments,  the  flatpack  gage  is  now  a  usable  tool  for  many  blast-load 
measurement  problems. 

Further  improvements  could  be  made.  First,  transverse  strain  coupling  into  foil  ^d-type 
elements  may  be  insignificanL  This  possibility  should  be  studied  further  (with  computations) 
because  it  would  allow  one  of  the  elements  to  be  eliminated.  Second,  the  mechanic^  and 
piezoresistance  properties  of  the  commercially  available  elements  should  be  measured.  Third,  the 
resfxjnse  time  (natural  period)  of  the  flatpack  gage  is  probably  adequate  for  most  explosive  load 
situations,  but  additional  design  changes  (such  as  mica  insulation  and  a  titanium  cover)  may 
improve  this  response  time.  Finally,  only  a  few  measurements  of  close-in  blast  pressures  have 
been  made  with  the  current  flatpack  gage  design,  and  although  the  results  are  good,  further  use  of 
the  gage  (and  comparisons  with  calciUations  and  other  measurements)  would  further  substantiate 
its  capabilities. 
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APPENDIX  A 

PIEZORESISTANCE  MATERIAL  ^  ELS 


The  material  models  for  ytterbium,  manganin,  and  constantan  used  in  this  project  are  based 
on  the  publications  of  Gupta  and  associates. This  Appendix  compares  model  inedictions  with 
test  data.  Table  A-1  lists  the  values  for  the  constitutive  and  piezoresistance  parameters  used  in 
PEEZOR.  Tensile  stresses  and  strains  are  positive-throughouL 

Figure  A-1  shows  stress-strain  curves  and  resistance  change  for  static  uniaxial  tension 

tests.  The  odd  nature  of  the  versus  ei  curve  exhibited  by  ytterbium  [Figure  A-l(b)l  arises 

because  of  ytterbium's  very  high  stress  sensitivity  and  its  low  strength.  During  the  elastic  phase  of 
the  tensile  test  the  resistance  change  is 


^  =  [(a  +  2p)  E  +  (1  -h  2v)l  ei  (A-1) 

where  the  first  term  represents  the  effect  of  the  axial  stress  and  the  second  term  accounts  for  the 
strains.  Using  the  values  for  a,  p,  E  and  v  in  Table  A-1,  the  first  term  in  the  brackets  has  a  value 

-5.05  and  the  second  term  1.73.  Thus,  the  slope  of  the  ^  versus  ej  curve  is  -3.32.  In  contrast, 

for  large  plastic  strains  the  stress  is  limited  to  about  1  kbar  [Figure  A- 1(a)]  and  the  ytterbium  is 
approximately  incompressible.  Thus,  ignoring  elastic  strains 


AR  1 

^  =  (a  +  2p)  (1  kbar)  +  (l  +2x^)Ei 


(A-2) 


Now  the  stress  term  has  a  value  of  about  -0.04  and  the  strain  coefficient  is  2.  The  same  effects 
occur  in  the  response  of  manganin  and  constantan,  but  for  these  materials  the  stress  sensitivity  is 


overwhelmed  by  the  strains  so  the  slopes  of  their 
strains. 


Ro 


versus  Ei  curves  approach  2  at  much  lower 
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Table  A-1.  Material  Properties 


E 

(kbar) 

[GPa] 

u 

Oy 

(kbar) 

[GPa] 

H 

(ld)ar) 

[GPa] 

•  * 

a 

(kbar"*) 

rGPa-'‘l 

P 

(kbari) 

[GPa-"*] 

A 

g 

Constantan 

1485 

[148.5] 

0.33 

8.2 

[0.82] 

409 

[40.9] 

-0.49  X  10-^ 
[-0.49  X  10-3] 

1.89  X  10-^ 

[1 .89  X  10-3] 

0 

0 

Manganin 

876 

[87.6] 

0.374 

2.3 

[0.23] 

24 

[2.4] 

-0.44  X  10-3 
[-0.44x10-2] 

-0.31  X  10-3 
[-0.31  X  10-2] 

0.0105 

60 

Ytterbium 

120 

[12.0] 

0.365 

0.49 

[0.049] 

212.5  (8.0  X  10-^) 
[21.25] 

122.2(1.7x10-^) 

-0.788  X  10-2 
[-0.788  X  10-”*] 

-1.71  X  10-2 
[-1.71  X  lO-”*] 

0.0425 

80 

[12.22] 

53.85  (3.0  X  10-^) 
[5.385] 

18.75  (5.4  X  10-3) 
[1.875] 

9.091  (7.6  X  10*3) 
[0.9091] 

0.8772  (1 .33  X  10-2) 
[0.08772] 


*For  ytterbium  a  multilinear  stress-strain  cunre  is  used.  The  quantities  in  parentheses  are  the 
values  of  plastic  strain  at  the  beginning  of  each  linear  segment. 

**For  ytterbium  the  pressure  dependence  of  a  is  -7.6  x  10-^  kbar2  [-7.6  x  10-2  GPa'2]. 
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Figure  A- 1 .  Static  uniaxial  tension  tests. 

Points  art  data  from  Rataranca  1.  Curvaa  ara  PIEZOR  modal  pradiction. 


Figure  A-2  shows  resistance  change  versus  nonna]  stress  for  uniaxial  strain  compression. 
The  PIEZOR  program  uses  the  ytterbium  curve  wi  Aout  Ac  pressure-dependent  coefficients  terms, 
so  it  should  be  modified  for  use  above  2  kbar  (200  MPa). 

Figure  A-3  shows  exponential  fits  to  Ae  piezoresistance  strain-hardenmg  coefficient  ("n)  for 
shock  wave  loading.  The  fits  fort)  could  be  improved  for  plastic  strains  less  Aan  1%. 


The  coefficient  in  equation  (11)  (Section  3)  is  defined  as  Ae  sensitivity  of  Ae  material's 

resistance  to  an  mcrement  of  plastic  strain.  The  determination  of  T)  m  References  3  and  4, 
however,  is  based  on  Ae  average  sensitivity  over  a  finite  step  m  plastic  strain.  Put  anoAer  way, 

AR  T5 

equation  (1 1)  calls  for  a  tangent  value  of  vs  curve,  whereas  References  3  and  4  give 

secant  values.  Tangent  values  were  determined  from  Ae  results  presented  m  References  3  and  4  as 
follows. 


First,  we  use  Ae  definitions  for  Ac  secant  values  of  Ae  ^  vs  curve  (denoted  ris)  and 

Ae  tangent  values  (denoted  tiO 

AR 

(A-3) 

yp 

K(Yp)dYp 

(A-4) 

Vs  _n 

where  Yp  =  "2"  Differentiating  (A-4)  and  substituting  from  (A-3)  gives 

If 

^[nsYp]=%Yp  +  % 

(A-5) 

Next,  Ae  vs  Yp  curve  is  constructed  fiom  Ae  results  given  m  References  3  and  4  usmg  (A~3), 
and  this  curve  is  fit  wiA  Ae  form 


f  =  A(i-c-«p) 

(A-6) 

For  Ais  form,  T|s  and  Tjt  are  given  by 

II 

l> 

1 

ft 

1 

(A-7) 

yp 

and 

Tit  =  A  g 

(A-8) 
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APPENDIX  B 

PIEZOR  COMPUTER  PROGRAM 


This  Appendix  contains  a  listing  of  the  PIEZOR  computer  progi^  and  its  associated 
subroutines,  lliese  are  followed  by  a  sample  input  set  and  a  paitid  listing  of  output 

The  &st  line  of  input  is  a  title.  Each  element  is  described  separately,  beginning  with  a 
material  name  (MAT  N).  Table  B-1  gives  the  definitions  of  the  input  parameters.  The  orientations 
of  the  elements  are  fixed;  elements  1  and  2  are  aligned  in  the  x  direction,  and  element  3  is  aligned  in 
the  y  direction. 
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Table  B-1.  PIEZOR  Variable  Deflnitions 


MAT 

Material  name 

MODELETA 

=  0  for  “n  =  constant 

=  1  for  11  =  Age* 

MODELH 

=  0  for  H  =  constant 

=  lforH(?) 

EM 

Elastic  modulus 

ANU 

Poisson's  ratio 

Yo 

Yield  stress 

H 

Hardening  modulus  (for  MODELH  =  0) 

NH 

EBARP,H 

Number  of  points  for  hardening  curve  (for  MODELH  =  1 

?,  H  pairs  for  hardening  curve 

ALPHA 

a 

BETA 

P 

ETAl 

11  (for  MODELETA  =  0)  or  A  (for  MODELETA  =  1 ) 

ETA2 

g  (for  MODELETA  =  1) 

LAMN  (i,i) 

II 

B-2 


PROGRAM  PIEZOR3 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


****  MODIFIED  6/14/91  TO  ADD  FEATURE  OF  H  AS  A  FUNCTION  OF  EPSBARP 
-  FEBRUARY,  1991 

PIEZOR3  CORRECTS  PROBLEM  WHERE  SIGMABAR  (SB)  DOES  NOT  TRACK  THE  YIELD 
ACCURATELY.  PIEZOR3  CALCULATES  B  MATRIX  AND  N, 

WHERE  N  -  (H+3*G*SBOLD/SBN) / (H+3*G) 

PROGRAM  ITERATES  UNTIL  SBN  CONVERGES. 

PIEZOR3  USES  SUBROUTINE  MATRXM2 

GIVEN  3  GAGE  HISTORIES  IN  DR/RO  VS  TIME,  WHERE  THE  GAGES  CONSIST 
OF: 

1.  A  STRESS  GAGE  IN  THE  RADIAL  DIRECTION 

2.  A  STRAIN  GAGE  IN  THE  RADIAL  DIRECTION 

3.  A  STRAIN  GAGE  PERPENDICULAR  TO  THE  STRESS  GAGE 

(IN  THE  CIRCUMFERENTIAL  DIRECTION) 

THIS  PROGRAM  CALCULATES  THE  NORMAL  STRESS  (SIGMAZ) ,  RADIAL  STRAIN, 

&  CIRCUMFERENTIAL  STRAIN  FOR  THESE  GAGE  DATA.  OTHER  QUANTITIES  ALSO 
CALCULATED  ARE  THE  REMAINING  STRESSES  AND  STRAINS  FOR  UP  TO  3 
MATERIALS  TO  COMPLETE  THE  SET  OF  STRESSES  £  STRAINS  IN  3  DIRECTIONS. 

UNITS:  OUTPUT  CONSISTENT  WITH  INPUT  UNITS 


INPUT : 

*  DATA  FILE  FOR  STRESS  GAGE 

*  DATA  FILES  (2)  FOR  STRAIN  GAGES 

*  MATERIAL  PROPERTIES  FOR  EACH  GAGE: 

AK  BULK  .MODULUS 

G  SHEAR  MODULUS 

YIELD  YIELD  STRESS 

H  WORK  HARDENING  MODULUS 

MODELH  INDIC-ATOR  FOR  H 

-  0  FOR  CONSTANT  H 

-  1  FOR  H  AS  A  FUNCTION  OF  EPSBARP 

INPUT  ARE  PAIRS  OF  POINTS  (EPSBARP, H). 

H  IS  CONSTANT  WITHIN  AN  EPSBARP  INTERVAL. 

*  PIEZORESISTIVE  CONSTANTS  FOR  EACH  GAGE 

ALPHA  PIEZORESISTIVE  COEFFICIENT  (1/Kbar) 

BETA  PIEZORESISTIVE  COEFFICIENT  (1/Kbar) 

MODELETA  INDICATOR  FOR  MODEL  TO  BE  USED  TO  CALCULATE  ETA,  THE 
PiiASTIC  DEFORMATION  COEFFICIENT 

-  0  FOR  ETA-ETAl  (A  CONSTANT) 

-  1  FOR  ETA-ETA1*ETA2*EXP(-ETA2*SQRT (0.75) *EBP) 

ETAl  CONSTANT  TERM  FOR  ETA  -  CONSTANT; 

ONE  OF  2  COEFFICIENTS  FOR  MODELETA-1 
ETA2  2ND  COEFFICIENT  FOR  MODELETA-1 

*  DIAGONAL  OF  LAMBDA,  THE  COUPLING  MATRIX  FOR  EACH  GAGE 

LAMBDA(1,1),  LAMBDA(2,2),  £  LAMBDA(3,3)  (0  <-  LAMBDA(3,3)  <-  1) 


OUTPUT : 

*  FILE  FOR006  FOR  BATCH  JOB 

FOR  EACH  GAGE  :  S,  EP,  SB,  SD,  DEPD,  EBP,  Y,  COND 

*  FILE  FOR008  -  BINARY  SCRATCH  FILE  FOR  STRESS-STRAIN  STATES 

*  FILE  FOR009  -  HISTORIES  FOR  STRESSES,  STRAINS,  EFFECTIVE  STRESS 

£  EFFECTIVE  PLASTIC  STRAIN  FOR  EACH  GAGE 
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PARAMETER  (MAXNPT-10001) 

IMPLICIT  REAL*8  (A-H,0-Z) 

INTEGER  COND,CONDS 
REAL*8  LAM,LAMN 

REAL*4  TS (MAXNPT) , RS (MAXNPT) , REl (MAXNPT) , RE2 (MAXNPT) , TEl (MAXNPT) 
1  ,TE2 (MAXNPT) ,DTS,DTE1,DTE2,TMIN,TMAX 
DIMENSION  AM(3,3,3) ,AMM(3,3) ,R(3) ,LAM(3, 3) , 

1  AMS (3, 3, 3), PS (3, 3) ,PM(1,3),QS(3,3),QM(1,3),QP{1,3) ,WS(3,3) , 

2  SS(3,3) ,SS2(3,1) ,TST(3,3) ,TS2{3,1) ,ROR(3) ,THETA(3) ,GM{3,3,3) , 

3  PROD (3, 3) ,FM(3) ,FM2 (3, 1) , OM (3, 3) ,B(3,3) 

DIMENSION  S(3,3) ,EP(3,3) ,SD(3,3),SB(3) ,EBP(3) ,DEBP(3) , 

1  Y(3) ,SOLD(3,3) ,P(3) ,EPOLD(3,3),SBOLD(3),EBPOLD(3) ,yOLD(3) , 

2  SDOLD(3,3) ,COND(3) , FR (3) , RT (3) ,DSD (3, 3) , ETA (3) , SDN (3, 3) , 

3  DEPD(3,3) ,SBN(3) ,SBNL(3) ,AN(3) 

CHARACTER* 80  TITLE 
CHARACTER* 50  FI LENAMEO , F I LENAMEH 
CHARACTER* 10  IDAT,MATNAM 
LOGICAL*!  LBATCH,LPC,LMATRX 

COMMON  IDAT, TITLE, lOUT, LBATCH, FILENAMEO, FILENAMEH 
C  LAMN(I,M)  -  LAMBDA  DIAGONAL  FOR  MATERIAL  M 

COMMON  /MAT/  AK (3) , G (3) , YO (3) , H (3)  , ETAl (3)  , ETA2 (3) , ALPHA (3) , 

1  BETA(3) ,MATNAM(3) ,LAMN(3,3) ,M0DELETA(3) ,M0DELH{3) ,NH(3) , 

2  EPH(50,3) ,HH(50,3) 

EQUIVALENCE  (FM ( 1) , FM2 (1,1)) 

DATA  MATNAM  /'MATERIAL  1',' MATERIAL  2', 'MATERIAL  3'/ 

DATA  DEGPRAD  /57. 295  (795/ 

DATA  THETA  /2*0,,90./ 

DATA  TOLER  /l.OE-5/ 


DEFINE  ETA  FOR  MODELETA-1 

ETAJ (Cl, C2, EBP) -C1*C2*EXP (-C2*SQRT (0.75) *EBP) 

DO  M-1,3 

THETA (M) -THETA (M) /DEGPRAD 
ENDDO 

READ  INPUT  DATA 

CALL  INPUT (MAXNPT, TS,RS,DTS,NPTS, TEl, REl, DTE1,NPTE1,TE2,RE2,DTE2, 
1  NPTE2,TMIN, TMAX) 

OPEN  SCRATCH  FILE 

OPiiN  (UNIT-8,  FORM-' UNFORMATTED'  ,  STATUS-' SCRATCH'  ) 

OPEN  (UNIT-8,  FORM-' UNFORMATTED' , STATUS-' NEW' ) 

INITIALIZE 

K-0 

KBEGl-1 
KBEG2-1 
NPT-0 
TIM-TMIN 
DTSNEW-0 . 

DO  M-1,3 
DO  1-1,3 
S(I,M)-0. 
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EP(I,M)-0. 

SD(I,M)-0. 

SOLD(I,M)-0. 

SDOLD(I,M)-0. 

EPOLD(I,M)-0. 

DSD(I,M)-0. 

SDN(I,M)-0. 

DEPD(I,M)-0. 

ENDDO 

SB(M)-0. 

SBOLD (M) -0 . 

EBP (M) -0 . 

DEBP (M) -0 . 

EBPOLD (M) -0 . 

SBN (M) >0 . 

SBNL (M) «0 . 

Y  (M)  -YO  (M) 

YOLD  (M)  -YO  (M) 

MATERIAL  IS  ASSUMED  ELASTIC  INITIALLY 

1  -  ELASTIC 

2  -  PLASTIC 

3  -  CURRENT  STEP  IS  ELASTIC,  NEXT  STEP  PLASTIC 

(FRACTION  OF  NORMAL  STEP  SIZE  FOR  TREATMENT  AT  CORNER) 

COND (M) -1 
ENDDO 


SET  UP  P*  &  Q* 

DO  M-1,3 

PS  (M,  1) -ALPHA  (M) 

PS(M,2)-1. 

PS<M,3)— 1. 

QS  (M,  1 )  -ALPHA  (M)  -f 2 .  *BETA  (M) 
QS  (M,  2) -ALPHA  (M) 

QS  (M,3)— 1. 

ENDDO 


SET  UP  ROTATION  MATRIX  OM  &  COUPLING  MATRIX  LAM 
DO  M-1,3 
DO  1-1,3 
DO  J-1,3 

LAM(I, J)-0. 

OM(I, J)-0. 

IF  (J  .EQ,  I)  THEN 

LAM ( I , J) -LAMN ( I , M) 

IF  (I  .EQ.  1)  THEN 
OM(I, J)-l. 

ELSE 

OM ( I , J) -COS ( THETA (M) ) *  *2 
ENDIF 
ENDIF 
ENDDO 
ENDDO 

OM  (2, 3) -SIN  (THETA  (M) )  **2 
OM  (3,  2) -SIN  (THETA  (M)  )  **2 
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CALL  MATMPy(LAM,3,3,OM,3,PROD) 

DO  1-1,3 
DO  J-1,3 

GM(I, J,M)-PROD(I, J) 

ENDDO 

ENDDO 

ENDDO 

C  LOOP  OVER  GAGE  DATA 

20  K-K+1 

IF  (K  .GT.  NPTS)  GO  TO  300 
IF  (TS(K)  .LT.  TMIN-1.E-4*DTS)  GO  TO  20 
IF  (TS(K)  .GT.  TMAX+1.E-4*DTS)  GO  TO  300 
IF  (NPT  .EQ.  0)  THEN 
R(1>-RS(K) 

GO  TO  25 
END  IF 

R(1)-RS(K)-RS(K-1) 

DTSNEW-DTS 

25  DO  I-KBEG1,NPTE1 

IF  (ABS(TE1{I)-TS(K))  .LE.  DTEl/2.)  THEN 
KBEGl-I 

IF  (NPT  .EQ.  0)  THEN 
R(2)-RE1 (I) 

ELSE 

R(2)-RE1<I)-RE1(I-1) 

ENDIF 
GO  TO  30 
ENDIF 
ENDDO 

30  DO  I-KBEG2,NPTE2 

IF  (ABS(TE2(I)-TS(K) )  .LE.  DTE2/2.)  THEN 
KBEG2-I 

IF  (NPT  .EQ.  0)  THEN 
R(3)-R£2 (I) 

ELSE 

R(3) -RE2 (I) -RE2 (I-l) 

ENDIF 
GO  TO  40 
ENDIF 
ENDDO 

40  CONTINUE 

C 

C  PREPARE  FOR  ITERATION,  DEPD(I,M)-0.  AS  1ST  GUESS 

42  ITER-0 

DO  M-1,3 
DO  1-1,3 

DEPD(I,M)-0. 

ENDDO 
SBNL (M) -0 . 

SBN(M)-0. 

ENDDO 

C 

C  ASSUME  MATERIAL  IS  ELASTIC  INITIALLY;  AN(M)-1  IF  ELASTIC 

DO  M-1,3 
COND (M) -1 
AN (M) -1.0 
FR(M)-0. 

RT (M) -0 . 
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ENDDO 

FRAC-l.OE+30 
NCORNER-0 
CONDS-0 
LMATRX-. FALSE. 

LPC-. FALSE. 

C 

C  CALCULATE  ETA  HERE 

45  CONTINUE 

DO  M-1,3 

IF  (MODELETA(M)  .EQ.  0)  THEN 
ETA (M) -ETAl (M) 

ELSE  IF  (MODELETA<M)  .EQ.  1)  THEN 

ETA (M) -ETAJ (ETAl (M) , ETA2 (M) , EBP (M) ) 

ELSE 

WRITE  (IOUT,9275)  M, MODELETA (M) 

9275  FORMAT ('  *****  ERROR  IN  SELECTION  OF  MODEL  FOR  ETA  — 

1  M  -'12,'  MODELETA  -'12/ 

2  '  ALLOWABLE  MODELS  ARE:'/ 

3  '  0  ETA-CONSTANT  OR' 

4  /'  1  ETA-ETA1*ETA2*EXP (-ETA2*EBP) ' ) 

STOP 

ENDIF 

ENDDO 

C 

C  CALCULATE  THE  M  4  M*  MATRICES  FOR  EACH  MATERIAL 
50  DO  M-1,3 

CALL  MATRXM2(M,AK(M),G(M),H(M),AN(M),AM) 

DO  J-1,2 

AMS(J,  1,M)-AM(J,3,M)  /AM(3,3,M) 

AMS(J,2,M)-AM(J,1,M)  -AM(J,3,M)  *AM{3,1,M)  /AM(3,3,M) 
AMS(J,3,M)-AM(J,2,M)  -AM(J,3,M)  *AM(3,2,M)  /AM(3,3,M) 

ENDDO 

AMS(3,1,M)-1./AM(3,3,M) 

AMS(3,2,M)— AM(3,1,M)  /AM{3,3,M) 

AMS(3,3,M)— AM{3,2,M)  /AM(3,3,M) 

ENDDO 

C 

C  CALCULATE  W*  4  B 

60  DO  M-1,3 

DO  1-1,3 

QM(1,I)-QS (M,I) 

PM(1,I)-PS(M,  I) 

DO  J-1,3 

AMM(I,  J)-AMS(I,  J,M) 

PRODd,  J)-GM(I,  J,M) 

ENDDO 

ENDDO 

CALL  MATMPY{PM, l,3,PROD,3,QP) 

CALL  MATMPY(QM, 1,3,AMM,3,PM) 

CALL  MATMPY(PM, l,3,PROD,3,QM) 

DO  J-1,3 

WS (M, J) -QP ( 1 , J) +QM ( 1 , J) 

ENDDO 


B-7 


on  non  on 


SS2 (1,1)- (AN (M) -1 . '  *SDOLD (3,M) 

SS2(2,l)-0. 

SS2(3,l)-0. 

CALL  MATMPY(AMM, 3,3,SS2,1,TS2) 

DO  1-1,3 

IF  (I  .NE.  3  THEN 

B (I,M) - (AN(M) -1 . ) *SDOLD {I,M) -TS2 (1,1) 

ELSE 

B(I,M)  — TS2(I,1) 

ENDIF 

ENDDO 

ENDDO 

CALCULATE  S* 

DO  M-1,3 
DO  1-1,3 

QM(1,1)-QS(M  I) 

SS2(I,1)-B(I,M) 

ENDDO 

CALL  MATMPY(QM,  1,3,SS2,1,TS2(M,  D) 

FM(M)-R(M)  -  SQRT(0.75)*ETA(M)*(EBP(M)-EBPOLD(M))  -  TS2(M, 1) 
ENDDO 

CALL  SIMQ(WS,FM, 3,IFLAG) 

CALCULATE  S  &  T  MATRICES 

DO  M-1,3 
DO  1-1,3 

DO  J-1,3 

PRODd,  J)-GM(1,  J,M) 

AMM(I,  J)-AMS(I,  J,M) 

ENDDO 

ENDDO 

CALL  MATMPY(PROD,3,3,FM2,l,SS2) 

CALL  MATMPY(AMM,3,3,SS2,1,TS2) 

DO  J-1,3 

SS(J,M)-SS2(J,  1) 

TST(J,M)-TS2(J,1)  +  B(J,M) 

ENDDO 

ENDDO 


COMPUTE  STRESS  DEVIATORS  &  STRAIN  DEVIATOR  INCREMENT 
DO  M-1,3 

S(l,M)-SOLD(l,M)  +TST(1,M) 

S(2,M)-SOLD(2,M)  +TST(2,M) 

S(3,M)-SOLD(3,M)  +SS(1,M) 

P(M)-(S(1,M)+S(2,M)  +S(3,M)  )/3. 

DO  1-1,3 

SD(1,M)-S(I,M)  -P  ;m) 

DSD(I,M)-SD(I,M)-SDOLD(I,M) 

ENDDO 

SB (M) -SORT (1.5* (SD(1,M) **2+SD(2,M) **2+SD(3,M) **2) ) 
ENDDO 
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TEST  FOR  YIELDING,  SET  MATERIAL  CONDITIONS 

IF  1ST  POINT  IS  PLASTIC  OR  A  CORNER,  STOP  CALCULATION 

-  MODIFIED  TEST  FOR  PLASTIC  CONDITION 

IF  (LMATRX)  GO  TO  95 

DO  70  M-1,3 
ELASTIC 

IF  (SB(M)  .LT.  YOLD (M) ) • THEN 
COND (M) -1 
PLASTIC 

ELSE  IF  (SB(M)  .GT,  YOLD (M)  .AND.  ABS (SBOLD (M) -YOLD (M) )  .LT. 

1  l.E-4*YOLD(M) )  THEN 

ELSE  IF  (SB(M)  .GT.  YOLD (M)  .AND.  SBOLD (M)  .GT.  YOLD (M) *0 . 9999) 
1  THEN 

COND (M) -2 

STRAIN  DEVIATOR  INCREMENT 

DEPD(1,M)-2.*SS(2,M) /3.-(SS(3,M) +TST(3,M) ) /3. 
DEPD{2,M)-2.*SS{3,M) /3 .- (TST (3,M) +SS{2,M) )/3. 

DEPD (3,M) -2 . *TST (3,M) /3 .- (SS {2,M) +SS (3,M) ) /3 . 

DO  1-1,3 

SON  ( I ,  M)  -SDOLD  ( I ,  M)  +2  .  *G  (M)  *DEPD  ( I ,  M) 

ENDDO 

SBN(M)-SQRT(1.5*(SDN(1,M)**2+SDN{2,M)**2+SDN(3,M)**2) ) 

AN  (M)  -  (H  (M)  +3 .  *G  (M)  *SBOLD  (M)  /SBN  (M)  )  /  (H  (M)  +3 .  *G  (M) ) 
CORNER 
ELSE 

COND  (M)  -3 
NCORNER-NCORNER+ 1 

LNDIF 

IF  (NPT  .EQ.  0  .AND.  COND (M)  .GT.  1)  GO  TO  950 
CONDS-CONDS+COND (M) 

CONTINUE 


IF  ONE  OR  MORE  MATERIAL  IS  PLASTIC  WHEN  ANOTHER  HAS  REACHED  A 
CJRNER,  RECALCULATE  M,  M*  AND  W  MATRICES  WITH  PLASTIC  CONDITION 
SO  THAT  THE  CORNER  CAN  BE  CALCULATED  MORE  ACCURATELY 


IF  (CONDS  .GE,  6  .AND.  CONDS  .LE.  8)  THEN 

IF  ({CONDS  .LE.  7  .AND.  NCORNER  .EQ.  1)  .OR. 

1  (CONDS  .EQ.  8  .AND.  NCORNER  .EQ.  2))  THEN 

IF  (.NOT.  LPC)  THEN 
LPC-.TRUE. 

CONDS-0 
NCORNER- 0 
GO  TO  50 

ENDIF 

ENDIF 

ENDIF 

C  IF  CORNER,  CALCULATE  FRACTION 

DO  M-1,3 

IF  (COND(M)  .EQ.  3)  THEN 
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Cl-1 .5*  (DSD  (1,M)  **2+DSD  (2,M)  **2+DSD  (3,M)  **2) 

C2-3* {SD0LD{1,M) *DSD(1,M)+SD0LD(2,M) *DSD(2,M) +SD0LD(3,M) 
*DSD(3,M) ) 

C3-SB0IiD  (M)  **2-Y0LD  (M) 

RAD-SQRT (C2**2-4 . *C1*C3) 

IF  (C3  .EQ.  0.)  THEN 
COND (M) -2 
CONDS-CONDS-1 
NCORNER-NCORNER- 1 
GO  TO  72 
ENDIF 

RT(l)-{-C2+RAD) / (2.*C1) 

RT(2)-{-C2-RAD) /{2.*C1) 

C  INSURE  ONE  OF  ROOT  MUST  BE  BETWEEN  0  AND  1 

DO  1-1,2 

IF  (ABS(RT(I))  ,LT.  l.OE-5)  RT(I)-0. 

IF  (ABS{RT(I)-1.0D0)  .LT.  l.OE-5)  RT(I)-1.0 
ENDDO 

IF  ((RT(1)  .LT.  0.  .OR.  RT{1)  .GT.  1.)  .AND.  (RT(2)  .LT. 

1  0.  .OR.  RT(2)  .GT.  1.))  THEN 

WRITE  (IOUT,9055)  M,COND (M) ,NPT, TIM,DTSNEW, RT (1) , RT (2) 
9055  FORMAT (//'  *****  STOP  —  ROOTS  OUTSIDE  RANGE  (0,1)  '/ 

1  '  M  -'12,'  COND  -'12,'  NPT  -'16,'  TIM,DTSNEW  -' 

2  1P2E12.4,'  RT(1),RT(2)  -'2E15.6) 

STOP  '  ***  ROOTS  OUTSIDE  RANGE  (0,1)' 

ENDIF 
DO  1-1,2 

IF  (RT(I)  .GT.  0.  .AND.  RT(I)  .LE.  1.)  THEN 
FR(M)-RT(I) 

GO  TO  72 
ENDIF 
ENDDO 

ENDIF 

72  CONTINUE 

ENDDO 

GO  TO  (960,960,95,75,80,80,80,80,80)  CONDS 
C  ELASTIC/PLASTIC 

75  LMATRX  -  .TRUE. 

GO  TO  50 

80  IF  (NCORNER  .EQ.  0)  GO  TO  75 

C  CORNER(S) 

DO  M-1,3 

IF  (COND(M)  .EQ.  3)  FRAC-MIN(FR(M) ,FRAC) 

ENDDO 

90  DTSNEW»DTSNEW*FRAC 

DO  J-1,3 

ROR(J)-R(J) 

R(J)-FRAC*R(J) 

ENDDO 

LMATRX  -  .TRUE, 

GO  TO  50 

C 

C  MAKE  CORRECTION  TO  STAY  ON  YIELD  SURFACE 

C 

95  CONTINUE 
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C  COMPUTE  STRAIN  DEVIATORS  AND  YIELD 
DO  M-1,3 

EP (l,M)-EPOLD(l,M) +SS (2,M) 

EP(2,M)-EPOLD(2,M) +SS(3,M) 

EP (3,M) -EPOLD(3,M)+TST(3,M) 

IF  (COND(M)  .EQ.  2)  THEN 

DEPD(l,M)-2.*SS(2,M)/3.-(SS(3,M) +TST(3,M) ) /3. 

DEPD (2,M) -2 . *SS (3> M) /3 (TST (3,M) +SS (2,M) ) /3 . 

DEPD(3,M) -2.*TST(3,M) /3.-(SS (2,M) +SS (3,M) ) /3. 

DO  1-1,3 

SDN ( I , M) -SDOLD ( I , M) +2 . *G (M) *DEPD ( I , M) 

ENDDO 

SBNL (M) -SEN (M) 

SEN (M) -SORT (1.5* (SDN (1,M) **2+SDN {2,M) **2+SDN(3,M) **2) ) 
AN  (M)  -  (H  (M)  +3  .  *G  (M)  *SEOLD  (M)  /SEN  (M)  )  /  (H  (M)  +3  .  *G  (M)  ) 

"  DEEP  (M)  -  (SEN  (M)  -SEOLD  (M)  )  /  (H  (M)  +3  .  *G  (M)  ) 

ELSE 

DEEP (M) -0 . 

END  IF 

EEP  (M)  -EEPOLD  (M)  -t-DEEP  (M) 

Y  (M)  -YOLD  (M)  +H  (M)  *DEEP  (M) 

ENDDO 

C 

C  TEST  FOR  CONVERGENCE  (SEN) 

C 

DO  M-1,3 

DELSEN-SEN (M) -SENL (M) 

IF  (AES(DELSEN)  .GT.  MAX (AES (TOLER*SEN (M) ) , 0 .DO) )  GO  TO  120 
ENDDO 
GO  TO  150 

120  MLAST-M 

IF  (ITER  .LE.  20)  THEN 
ITER-ITER+1 
DO  M-1,MLAST 

SENL (M) -SEN (M) 

ENDDO 
GO  TO  45 
ELSE 

WRITE  (IOUT,9050)  NPT, K, TIM, TOLER, (M, ETA (M) , 

1  SENL (M) , SEN (M) ,  M-1,MLAST) 

9050  FORMAT('  *****  STOP:  ITER>20,  NPT,K-'2I5,'  TIM  -' IPEll . 3, 

1  '  TOLER  -'E10.3/('  M  -'13,'  ETA  -'IPEll. 3, 

2  '  SENL,  SEN  -'2E13.5)) 

STOP  '  ITERATION  EXCEEDED  20' 

END  IF 

C 

C 

C  UPDATE  STRESS  AND  STRAIN  STATES 

150  CONTINUE 
NPT-NPT+1 
TIM-TIM+DTSNEW 
DO  M-1,3 
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c 

c 

c 

200 

C 

C 

C 


c 

c 

300 

320 

1110 

1115 

340 

C 

C 

1120 


DO  J-1,3 

SOLD{J,M)-S(J,M) 

SDOLD (J,M)-SD(J,M) 

EPOLD ( J, M) -EP ( J, M) 

ENDDO 

SBOLD  (M)  -SB  (M) 

EBPOLD (M) -EBP (M) 

IF  {COND(M)  .EQ.  2)  YOLD (M) -SB (M) 

CALCULATE  H 

CALL  STRNHARD(M,NH(M)  ,EPH{1,M) ,HH(1,M)  ,MODELH (M) , EBP (M) , H (M) ) 
ENDDO 


WRITE  STRESSES  &  STRAINS  TO  BINARY  FILE  FOR008 

WRITE  (8)  NPT, TIM, ((S(I,M), 1-1,3), (SD(I,M), 1-1,3), SB(M),Y(M), 

1  P (M) , (SS (I,M) , 1-1,3) , (TST (I,M) , 1-1, 3) , (EP (I,M) , I-l, 3) , 

2  (DEPD(I,M),l-l,3),EBP(M),COND(M),  M-1,3) 

CHECK  TO  SEE  IF  TIME  STEP  ADJUSTMENT  HAS  BEEN  MADE 

IF  (NCORNER  .EQ.  0)  THEN 
GO  TO  NEXT  TIME  STEP 
GO  TO  20 
ELSE 

RESET  TIME  STEP  AND  DR/RO 
DTSNEW-DTS* (1 . -FRAC) 

DO  J-1,3 

R ( J) -ROR ( J) * (1 . -FRAC) 

ENDDO 
GO  TO  42 
ENDIF 


WRITE  OUTPUT  TO  FOR006  (BATCH) 

CONTINUE 
DO  340  MT-1,3 
REWIND  (8) 

READ  (8, END-340)  NPT, TIM, ( (S (I,M) , 1-1, 3) , (SD (I,M) , I-l, 3) , SB (M) , 

1  Y(M),P(M), (SS{I,M) ,1-1,3) , (TST(1,M),I-1,3), (EP (I,M) , I-l, 3) , 

2  (DEPD(I,M),I-l,3),EBP(M),COND(M),  M-1,3) 

IF  (NPT  .EQ.  1)  WRITE  (IOUT,1110)  TITLE,  MT, MATNAM (MT) 
FORMATdHl, /IX,  A79//30X, 'MATERIAL  ',11,'  -  'AlO 

1  //'  I  TIME  St31  EPS  II]  EPS  12]' 

2  '  S[l]  S[2]  EPSI3]  P' 

3  '  SBAR  YIELD  C  EPSBARP'/) 

WRITE  (IOUT,1115)  NPT,TIM,S(3,MT),EP(1,MT),EP(2,MT),S(1,MT), 

1  S (2, MT) , EP (3, MT) , P (MT) , SB (MT) , Y (MT) ,  COND (MT) ,  EBP (MT) 
FORMAT(I6,1P8E11.3,2E10.3,I2,E10.3) 

GO  TO  320 
CONTINUE 


WRITE  STRESS-STRAIN  HISTORIES  TO  FOR009  FOR  POST  PROCESSING 

OPEN  (UNIT-9,  STATUS-' NEW' ,RECL-300) 

WRITE  (9,1120)  IDAT, TITLE 
FORMAT ( IX, A1 0 , A6 9 ) 

WRITE  (9,1130)  (M,MATNAM(M) ,  M-1,3) 
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1130  FORMAT(3(I5,'-',A10)) 

WRITE  (9  1135) 

1135  FORMAT (IX, 'TIME,  (SI (M) ,  S2 (M)  ,S3 (M)  ,EPS1 (M) ,EPS2 (M) ,EPS3 (M) , ' , 

1  ' EBP (M) , SEAR (M) ,  M-1,3)') 

REWIND  (8) 

350  READ  (8, END-370)  NPT, TIM, ( (S (I,M) , I-l, 3) , (SD (I,M) , I-l, 3) , SB (M) , 

1  y(M),P(M), (SS(I,M),I-1,3), (TST(I,M) ,1-1,3) , (EP (I,M) , I-l, 3) , 

2  (DEPD(I,M),I-l,3),EBP(M),COND(M),  M-1,3) 

WRITE  (9,1125)  TIM,  ( (S ( J, M) , J-1, 3) , (EP ( J,M) , J-1, 3) , EBP (M) , SB (M) , 

1  M-1,3) 

1125  FORMAT (1P25E12. 4) 

GO  TO  350 
370  CONTINUE 
CLOSE  (9) 

CLOSE  (8) 

CLOSE  (lOUT) 

GO  TO  990 

910  WRITE  (IOUT,9020)  FILENAMEH 

9020  FORMAT('  -***  ERROR  IN  OPENING  OUTPUT  FILE  ',A50) 

GO  TO  990 

950  WRITE  (IOUT,9057)  NPT,M,COND (M) 

9057  FORMAT ('  ****  ERROR  EXIT  —  NPT  -'12,'  MATERIAL  -'12,'  COND  -'12) 
GO  TO  990 

960  WRITE  (IOUT,9030)  CONDS 

9030  FORMAT ('  ****  .  ERROR  EXIT  —  COND (1) +COND (2) +COND (3)  -'12) 

990  STOP 
END 
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SUBROUTINE  INPUT (MAXNPT, TS , RS, DTS, NPT, TEl, REl , DTEl , NPTEl , TE2 , RE2 , 

1  DTE2,NPTE2,TMIN,TMAX) 

C 

C  READS  STRESS  VS  TIME,  STRAINS  VS  TIME,  &  MATERIAL  PROPERTIES 
C 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*4  TS (MAXNPT) , RS (MAXNPT) , REl (MAXNPT)  ,  RE2 (MAXNPT) , TEl (MAXNPT) , 

1  TE2 (MAXNPT) , DTS, DTEl, DTE2,TSH1FT,VSHIFT,T0A,TMIN,TMAX 
REALMS  LAMN 
DIMENSION  E(2),ANU(2) 

CHARACTER* 50  FI LENAME , F I LENAMEO , F I LENAMEH 
CHARACTER*80  TITLES, TITLEl, TITLE 
CHARACTER* 40  NAME 

CHARACTER* 10  IDAT , A1 , A2 , A3 , A4 , MATNAM 
LOGICAL*!  LBATCH 

COMMON  IDAT, TITLE, lOUT, LBATCH, FI LENAMEO, FILENAMEH 
C  LAMN(I,M)  -  LAMBDA  DIAGONAL  FOR  MATERIAL  M 

COMMON  /MAT/  AK (3) , G (3) , YO (3) , H (3)  ,ETA1 (3) , ETA2 (3) , ALPHA (3) , 

1  BETA(3) ,MATNAM(3) ,LAMN(3,3)  ,MODELETA(3) ,MODELH (3) , NH (3) , 

2  EPH{50,3) ,HH(50,3) 

C 

C  FOR  BATCH  JOBS,  MATERIAL  PROPERTIES  AND  FILENAMES  OF  GAGE  DATA 
C  ARE  READ  FROM  FOR005,  GAGE  DATA  FROM  FOROOl 

C 

CALL  DATE (IDAT) 

IN-S 

IOUT-6 

C 

C  WRITE  DATE 

WRITE  (IOUT,1100)  IDAT 

1100  FORMAT (IHl,'  DATE  -'AlO, 30X, 'BATCH  JOB') 

C 

C  READ  MATERIAL  PROPERTIES 

C  ***  MATERIAL  1  -  MATERIAL  OF  STRESS  GAGE 

C  ***  MATERIAL  2  -  MATERIAL  OF  STRAIN  GAGE  1 

C  ***  MATERIAL  3  -  MATERIAL  OF  STRAIN  GAGE  2 

READ  (IN, 1025)  TITLE 
1025  FORMAT (A80) 

WRITE  (IOUT,1110)  TITLE 
1110  FORMAT ( IX, A7 9) 

C 

C  MODELH  -  0  FOR  CONSTANT  H 

C  -  1  FOR  H  AS  A  FUNCTION  OF  EPSBARP 

C  (input  as  points;  H  is  constant  in  an  interval  of  EPSBARP) 

C 

DO  MAT-1,3 

READ  (IN, 1028)  Al, MATNAM (MAT ) ,A2,MODELETA (MAT) , A3, MODELH (MAT) 

1028  FORMAT (2A10, 2 (AlO, 110) ) 

WRITE  (IOUT,1028)  Al, MATNAM (MAT ) ,A2,MODELETA (MAT) , A3, 

1  MODELH (MAT) 

READ  (IN,1030)  A1,E{MAT),A2,ANU(MAT),A3,Y0(MAT),A4,H(MAT) 

WRITE  (lOUT,1035)  Al,  E  (MAT)  ,  A2,  ANU  (MAT)  ,  A3,  YO  (MAT)  ,  A4,  H  (MAT) 

C  READ  (IN,  1030)  Al,  AK  (MAT)  ,  A2,G  (MAT) ,  A3,  YO  (MAT)  ,  A4,  H  (MAT) 

1030  FORMAT(4(A10,E10.3)  ) 

C  WRITE  (IOUT,1035)  Al, AK (MAT) , A2,G (MAT) , A3, YO (MAT)  ,  A4,  H (MAT) 

IF  (MODELH (MAT)  .EQ.  1)  THEN 
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READ  (IN, 1033)  A1,NH{MAT) 

FORMAT (AlO, 110) 

WRITE  (IOUT,1033)  A1,NH(MAT) 

READ  (IN,1050)  Al, (EPH {I,MAT) ,HH (I,MAT) ,  I-1,NH(MAT)) 

FORMAT (AlO,  6E10.3/ (8E10.3) ) 

WRITE  (IOUT,1055)  Al, (EPH(I,MAT) ,HH(I,MAT) ,  I=1,NH(MAT)) 
FORMAT (AlO, 1P6E10 . 3/ (1P8E10 .3) ) 

END  IF 

FORMAT(4(AiC,lPE10.3) ) 

INITIALIZE  H 

IF  (MODELH(MAT)  ,EQ.  1) ■H(MAT)-HH(1,MAT) 

READ  (IN, 1030)  A1,ALPHA(MAT) ,A2,BETA(MAT) ,A3,ETA1 (MAT)  , 

A4,ETA2 (MAT) 

WRITE  (IOUT,1035)  Al, ALPHA (MAT) ,A2, BETA (MAT) , A3, ETAl (MAT) , 
A4,ETA2 (MAT) 

AK(MAT)-E (MAT) / (3 . * (1 . -2*ANU (MAT) ) ) 

G (MAT) -E (MAT) / (2 . * ( 1 . +ANU (MAT) ) ) 

WRITE  (IOUT,1040)  AK (MAT) , G (MAT) 

FORMAT ( 1 OX, 'AK  -'1PE12.4,'  G  -'E12.4) 

READ  (IN, 1030)  A1,LAMN(1,MAT) ,A2,LAMN(2,MAT) ,A3,LAMN(3,MAT) 
WRITE  (IOUT,1035)  A1,LAMN(1,MAT) ,A2,LAMN(2,MAT) ,A3,LAMN(3,MAT) 
ENDDO 

C 

C  BATCH  INPUT  USES  FREE-FIELD  FORMAT  FOR  READING  GAGE  DATA 
C 

CALL  READB (IN, lOUT, FILENAME, TITLES, NPT, DTS, TS, RS, TSHIFT, TOA) 

WRITE  (IOUT,1110)  FILENAME 

CALL  READB ( IN, lOUT , FILENAME, TITLEl , NPTEl , DTEl , TEl , REl , TSHIFT , TOA) 
IF  (DTEl  .NE.  DTS)  THEN 

WRITE  (IOUT,9900)  DTEl, DTS 

9900  FORMAT('  *****  DTEl  &  DTS  ARE  NOT  THE  SAME... DTEl  -'1PE11.3, 

1  '  DTS  -'Ell. 3/'  >»»  STOP  IN  SUBROUTINE  INPUT') 

STOP 

ENDIF 

WRITE  (IOUT,1110)  FILENAME 

CALL  READB ( IN, lOUT, FILENAME, TITLEl, NPTE2, DTE2, TE2 , RE2, TSHIFT, TOA) 
IF  (DTE2  .NE.  DTS)  THEN 

WRITE  (IOUT,9905)  DTE2,DTS 

9905  FORMAT ('  *****  DTE2  &  DTS  ARE  NOT  THE  SAME...DTE2  -'1PE11.3, 

1  '  DTS  -'El  1.3/'  »»>  STOP  IN  SUBROUTINE  INPUT') 

STOP 

ENDIF 

WRITE  (IOUT,1110)  FILENAME 
C 

C  SETUP  DATA 

C 

200  K-0 

TMIN-MAX(TS (1)  ,TE1 (1) ) 

TMIN-MAX (TMIN, TE2 (1) ) 

IF  (ABS(TMIN-TS(1) )  .LE.  l.E-4*DTS)  THEN 
TMIN-TS (1) 

GO  TO  210 

ELSE  IF  (ABS(TMIN-TE1(1))  .LE.  l.E-4*DTEl)  THEN 
TMIN-TEl (1) 

GO  TO  210 

ELSE  IF  (ABS (TMIN-TE2 (1) )  .LE.  l.E-4*DTE2)  THEN 
TMIN-TE2 (1) 


1033 

1050 

1055 

1035 

C 

1 

1 

1040 
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END  IF 

210  TMAX-MIN(TS(NPT),TE1(NPTE1)) 

TMAX-MIN (TMAX, TE2 (NPTE2) ) 

DO  I-1,NPTE1 

IF  (ABS (TEl (I) -TMAX)  .LE.  DTEl/2.)  GO  TO  220 
ENDDO 

220  NPTEl-I 

DO  I-1,NPTE2 

IF  (ABS(TE2(I)-TMAX)  .LE.  DTE2/2.)  GO  TO  230 
ENDDO 

230  NPTE2-I 
RETURN 

900  WRITE  (IOUT,9910)  FILENAME 

9910  FORMAT ('  *****  ERROR  IN  OPENING  FOROOl,  FILE  'A50) 
GO  TO  990 

910  WRITE  (IOUT,9920)  FILENAMEO 

9920  FORMAT ('  *****  ERROR  IN  OPENING  FILE  :'A50) 

990  STOP 
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SUBROUTINE  READB ( IN, lOUT, FILENAME, TITLE, NPT, DT, T, A, TSHIFT, TOA) 


READ  DATA  FROM  FOROOl  WITH  FREE-FIELD  FORMAT 
WRITE  ERROR  MESSAGE  TO  UNIT  ' lOUT' 

IMPLICIT  REAL*4  (A-H,0-Z) 

CHARACTER  FILENAME*50,  TITLE*80 
DIMENSION  T(1),A(1) 

C 

READ  (IN, 1010)  FILENAME 
1010  FORMAT (A50) 

OPEN  (UNIT-1,FILE-FILENAME, STATUS-' OLD' ,ERR-900) 
CALL  RDFLATPK (1 , TITLE, NPT, DT,  T,  A,  TSHIFT, TOA) 

CLOSE  (UNIT-1) 

RETURN 

900  WRITE  (IOUT,1100)  FILENAME 

1100  FORMAT ('  *****  ERROR  IN  OPENING  FILE  :'A50 

1  /'  »»>  STOP  IN  SUBROUTINE  READB') 

STOP 
END 
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SUBROUTINE  RDFLATPK (IN, TITLE, NPT,DT, T, A,TSHIFT, TOA) 


1030 


READ  DATA  WRITTEN  IN  FREE-FIELD  FORMAT  FROM  FILE  'IN' 

1 .  title 

2.  npt,  dt,  tshift,  toa;  (DATA  WILL  BE  SHIFTED  BY  T(l)) 

3.  (a(i),  i-l,npt)  (FREE-FIELD) 

IMPLICIT  REAL*4  (A-H,0-Z) 

DIMENSION  T(l) ,A(1) 

CHARACTER* 80  TITLE 
TSHIFT-0 . 

TOA-0 . 

READ  (IN, 1030)  TITLE 
FORMAT (ABO) 

READ  (IN,*)  NPT,  DT,  TSHIFT,  TOA 
READ  (IN,*)  (A(I),I-1,NPT) 

DO  I«1,NPT 

T(I)-(I-1) *DT  +  TSHIFT 
ENDDO 
RETURN 
END 
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c 

c 

c 

c 

c 

c 
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SUBROUTINE  MATMPY (A, N,M, B, L, C) 

MULTIPLIES  MATRICES  A  BY  B  AND  STORES  IN  C 
A  IS  n  X  m 
B  IS  m  X  1 
C  IS  n  X  1 

IMPLICIT  REAL* 8  (A-H,0-Z) 

DIMENSION  A(N, 1),B(M, 1),C(N, 1) 

DO  10  I-1,N 
DO  10  J-1,L 
C(I, J)-0. 

DO  10  K“1,M 

C(I,  J)-C(l,  J)+A(I,K)  *B  (K,  J) 

RETURN 

END 
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SUBROUTINE  MATRXM2  (M,AK,G,H,AN,AM) 

CALCULATE  THE  ELEMENTS  OF  M  —  FOR  PIE20R3 
NEW  DEFINITION  OF  M  MATRIX 

IMPLICIT  REAL*8  (A-H,0-2) 

DIMENSION  AMO,  3,1) 

INTEGER  COND 

AKP-AK+4 . *G*AN/3 . 

AKM-AK-2 . *G*AN/3 . 

DO  1-1,3 
DO  J-1,3 

IF  (J  .EQ.  I)  THEN 

AM(I,J,M)  -  AKP 

ELSE 

AM(I,J,M)  -  AKM 

END  IF 
ENDDO 
ENDDO 
RETURN 
END 
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SUBROUTINE  SIMQ(A,B,N,KS) 

IMPLICIT  REAL*8  {A-H,0-Z) 

DIMENSION  A{1),B(1) 

C  FORWARD  SOLUTION 

TOL-0 . 0 
KS-0 
JJ— N 

DO  65  J-1,N 
JY-J+1 
JJ-JJ+N+1 
BIGA-0 . 

IT=JJ-J 
DO  30  I-J,N 

C  SEARCH  FOR  MAXIMUM  COEFFICIENT  IN  COLUMN 

IJ-IT+I 

IF  (ABS(BIGA)-ABS(A(IJ) ) )  20,30,30 
20  BIGA-A(IJ) 

IMAX-I 
30  CONTINUE 

C  TEST  FOR  PIVOT  LESS  THAN  TOLERANCE  (SINGULAR  MATRIX) 

IF  (ABS (BIGA) -TOL)  35,35,40 
35  KS-1 
RETURN 

C  INTERCHANGE  ROWS  IF  NECESSARY 

40  Il-J+N*{J-2) 

IT-IMAX-J 
DO  50  K-J,N 

11- Il+N 

12- I1+IT 
SAVE-A(Il) 

A(I1)-A(I2) 

A (12) -SAVE 

C  DIVIDE  EQUATION  BY  LEADING  COEFFICIENT 

50  A(I1)-A(I1) /BIGA 
SAVE-B(IMAX) 

B(IMAX)-B(J) 

B(J)-SAVE/BIGA 

C  ELIMINATE  NEXT  VARIABLE 

IF  (J-N)  55,70,55 
55  IQS-N*(J-1) 

DO  65  IX-JY,N 

IXJ-IQS+IX 

IT-J-IX 

DO  60  JX-JY,N 

IXJX-N* (JX-1)+IX 

JJX-IXJX+IT 

60  A(IXJX)-a(IXJX)-(A(IXJ) *A(JJX) ) 

65  B(IX)-B(IX)-(B(J)*A(IXJ) ) 

C  BACK  SOLUTION 

70  NY-N-1 
IT-N*N 

DO  80  J-1,NY 
lA-IT-J 
IB-N-J 
IC-N 

DO  80  K-1,J 

B(IB)-B(IB)-A(IA)  ‘BdO 
lA-IA-N 
80  IC-IC-1 
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RETURN 

END 
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SUBROUTINE  STRNHARD (M,NH,EPH,HH,MODELH,EBP,H) 


C 

C  CALCULATE  H  BASED  ON  MODELH 
C 

C  MODELH  -  0  FOR  H  =  CONSTANT  (INPUT) 

C  1  FOR  H  AS  A  FUNCTION  OF  EPSBARP  (TABLE  OF  INPUT  VALUES) 

C  (H  IS  CONSTANT  IN  EACH  EPSBARP  INTERVAL) 

C 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  EPH(50) ,HH(50) 

IF  (MODELH  .EQ.  0)  RETURN 
IF  (MODELH  .EQ.  1)  THEN 
DO  I-1,NH 

IF  (EBP  .LE.  0.)  RETURN 
IF  (EBP  .LE.  EPH(I))  THEN 
H  -  HH(I) 

GO  TO  100 

ENDIF 

ENDDO 

IF  (EBP  .GT.  EPH(NH))  H-HH (NH) 

ENDIF 

100  RETURN 
END 
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SAMPLE  INPUT  DATA 


PLAS  CELL  -  3A  -  ETA{M)  NOT  0 , STRNHARDENING  (NIKE  DATA  OF  6-91)  ***PIEZOR3 
MAT  1  -  YTTERBIUM  MODELETA-  0  MODELH  -  1 

EM  -  120.0  ANU  -  0.365  YOY  -  0.49  HY  -  0. 

NH  -  6 

EBARP,H  -  8.000E-04  2.125E+02  1.700E-03  1.2222E+2  3.000E-03  5.3846E+1 
5.400E-03  1.875E+01  7.600E-03  9.0909E+0  1.330E-02  8.7719E-1 

ALPHAY  «  -7.88E-03  BETAY  -  -1.71E-02  ETAIY  -  1.13  ETA2Y  -  0. 

LAMN(1,1)  1.  LAMN(2,2)  1.  LAMN(3,3)  1. 

MAT  2  -  CONSTANTAN  MODELETA-  0 

EM  -  1484.4181  ANU  -  0.33  YOC  -  8.2  HC  -  409. 

ALPHAC  -  -0.49E-04  BETAC  -  1.89E-04  ETAIC  -  0.0 

LAMN(1,1)  1.  LAMN(2,2)  1.  LAMN(3,3)  1. 

MAT  3  -  CONSTANTAN  MODELETA-  0 

EM  -  1484.4181  ANU  -  0.33  YOC  -  8.2  HC  -  409. 

ALPHAC  -  -0.49E-04  BETAC  -  1.89E-04  ETAIC  -  0.0 

LAMN(1,1)  1.  LAMN(2,2)  1.  LAMN(3,3)  1. 

YBPZl.RES 

ClPZl.RES 

C2PZ1.RES 


YTTERBIUM  CELL  FOR  PIEZOR  TEST  DRl/RO 

51  2.0000E-02  O.OOOOE+00  O.OOOOE-t-00 

O.OOOOE+00  2.0664E-02  4.1323E-02  6.5196E-02  9.2118E-02  1.2103E-01 

1.5244E-01  1.8552E-01  2.1967E-01  2.5490E-01  2.9013E-01  2.9228E-01 

2.9487E-01  2.9776E-01  3.0088E-01  3.04i5E-01  3.0752E-01  3.1096E-01 

3.1444E-01  3.1795E-01  3.2147E-01  3.2609E-01  3.3071E-01  3.3533E-01 

3.3996E-01  3.4460E-01  3.4924E-01  3.5388E-01  3.5854E-01  3.6319E-01 

3.6785E-01  3.7252E-01  3.7719E-01  3.8187E-01  3.8656E-01  3.8986E-01 

3.9066E-01  3.9128E-01  3.9173E-01  3.9205E-01  3.9226E-01  3.8755E-01 

3.8284E-01  3.7813E-01  3.7344E-01  3.6874E-01  3.6405E-01  3.5937E-01 

3.5469E-01  3.5002E-01  3.4536E-01 


CONSTANTAN  CELL  1  FOR  PIEZOR  TEST  DRl/RO 

51  2.0000E-02  O.OOOOE+OO  O.OOOOE+00 


O.OOOOE+00  1.8270E-04  3.6546E-04 

1.0959E-03  1.2784E-03  1.4607E-03 

5.8483E-03  7.8579E-03  9.8663E-03 

1.7828E-02  1.9810E-02  2.1788E-02 

1.3781E-02  1.1776E-02  9.7702E-03 

1.7336E-03  -2.7880E-04  -2.2923E-03 

-1.0117E-02  -1.2042E-02  -1.3976E-02 

-1.3821E-02  -1.1799E-02  -9.7779E-03 

-1.7074E-03  3.0706E-04  2.3202E-03 


5.4803E-04  7.3067E-04  9.1319E-04 

1.6430E-03  1.8254E-03  3.8374E-03 

1.1864E-02  1.3856E-02  1.5844E-02 

1.9788E-02  1.7787E-02  1.5785E-02 

7.7629E-03  5.7544E-03  3.7446E-03 

-4.3070E-03  -6.2970E-03  -8.2017E-03 

-1.5919E-02  -1.7869E-02  -1.5844E-02 

-7.7582E-03  -5.7400E-03  -3.7232E-03 


CONSTANTAN  CELL  2  FOR  PIEZOR  TEST  DRl/RO 

51  2.0000E-02  O.OOOOE+OO  O.OOOOE+OO 


O.OOOOE+OO  l,8270E-04  3.6546E-04 

1.0959E-03  1.2784E-03  1.4607E-03 

1.0048E-03  5.9364E-04  1.8152E-04 

-3.1034E-04  -4.2254E-04  -5.3608E-04 

1.1318E-03  1.5467E-03  1.9609E-03 

3.6091E-03  4.0189E-03  4.4282E-03 

5.3434E-03  5.4224E-03  5.4995E-03 

4.8449E-03  4.4418E-03  4.0381E-03 

2.4131E-03  2.0049E-03  1.5958E-03 


5.4803E-04  7.3067E-04  9.1319E-04 

1.6430E-03  1.8254E-03  1.4155E-03 

2.1342E-05  -8.8538E-05  -1.9901E-04 

-1.1747E-04  2.9986E-04  7.1637E-04 

2.3739E-03  2.7865E-03  3.1982E-03 

4.8367E-03  5.1757E-03  5.2616E-03 

5.5750E-03  5.6488E-03  5.2474E-03 

3.6327E-03  3.2272E-03  2.8206E-03 
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APPENDIX  C 


EXPLOSIVE  BLAST  CALCULATIONS 


The  L2D  two-dimensional  finite-difference  code*  was  used  to  perform  two  calculations  to  help 
in  understanding  our  preliminary  experiments,  Tests  2  and  3  (see  Appendix  D)  and  to  help  design  later 
experiments.  T^e  mesh  layout  for  these  calculations  is  shown  in  Figure  C-1.  The  steel  was  modeled 
as  elastic  perfectly  plastic  with  a  flow  stress  of  10  kbar  (I  GPa).  A  planar  interface  was  prescribed 
within  the  steel  target  at  the  1.27-cm  depth,  corresponding  to  the  location  of  the  epoxy  layer  in  Test  2. 
In  the  first  calculation,  no  slip  was  allowed  at  the  interface;  in  the  second  calculation,  the  interface  was 
allowed  to  slip  without  friction.  These  conditions  bound  the  in-plane  properties  of  the  epoxy  layer, 
however,  the  effect  of  the  finite  thickness  of  the  epoxy  on  stresses  and  displacements  in  the  axid 
direction  were  not  modeled  in  the  calculations.  The  no-slip  interface  calculation  also  represents  the 
geometry  of  Test  3  with  the  omission  of  the  steel  and  epoxy  strips  comprising  the  flatpack  gages. 

Figures  C-2  through  C-4  show  the  normal  stresses  and  in-plane  strains  calculated  at  the  surface 
and  at  the  interface  for  the  no-slip  condition.  Except  for  the  peaks,  the  calculated  surface  and  interface 
stresses  at  r  =  2.54  cm  are  very  much  alike,  suggesting  that  a  gage  placed  at  a  non-slipping  interface 
1.27-cm  below  the  loaded  surface  should  provide  a  reascmably  accurate  measurement  of  the  surface 
load.  These  calculated  stresses  also  resemble  the  manganin  element  records  from  those  tests  in  that  the 

pulse  duration  is  about  25  p.s  and  the  peak  stress  at  the  surface  is  greater  than  at  the  interface. 

However,  both  calculated  peaks  are  considerably  lower  than  the  (approximate)  measurements. 

The  calculated  stresses  at  r  =  5.08  cm  are  also  similar  to  each  other,  and  the  ~1.5  kbar 
magnitude  is  consistent  with  the  interpretation  of  the  manganin  and  strain-gage  records. 

The  strains  from  the  no-slip  calculation,  shown  in  Figures  C-3  and  C-4,  are  very  different  from 
the  strain  gage  records.  First,  the  overall  levels  of  strain  are  about  an  order  of  magnitude  lower  than 
those  indicated  by  the  strain  gages.  Second,  the  large  oscillations  in  the  measurements  are  not 
predicted  in  the  calculations. 

Figure  C-5  shows  that  allowing  frictionless  sliding  at  the  interface  has  little  effect  on  the 
predicted  stresses  at  the  interface.  Figures  C-6  and  C-1  show  that  the  sliding  interface  condition  does 
change  the  sense  of  the  strains  at  the  interface  (strains  in  the  cylinder  at  the  interface  boundary),  but  the 
strain  levels  are  unchanged  and  bear  little  resemblance  to  the  strain-gage  records. 

The  probable  cause  for  the  general  lack  of  resemblance  between  the  calculated  and  measured 
strains  is  that  the  finite  thickness  of  the  epoxy  layer  was  not  included  in  the  calculation. 


*T.  Cooper,  "Users  Manual  for  SWE2D,"  Swedish  Detonic  Research  Foundation,  Report  DS  1982:14  (1982). 
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Figure  C-2.  Calculated  normal  stresses  with  no  slip  at  the  interface. 
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Figure  C-3.  Calculated  radial  strains  with  no  slip  at  the  interface. 
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Figure  C-4.  Calculated  circumferentiai  strains  with  no  slip  at  the  interface. 
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Figure  C-5.  Calculated  normal  stresses  at  the  interface  (1 .27-cm  depth). 
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Figure  C-6.  Calculated  racfial  strains  at  the  intertaoe  (1 .27-cm  depth). 
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Figure  C-7.  Calculated  circumferential  strains  at  the  interlace  (1 .27-cm  depth). 
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APPENDIX  D 


PRELIMINARY  EXPERIMENTS 


The  design  of  the  flatpack  evolved  during  the  project,  as  directed  by  developmental 
experiments  and  design  calculations.  The  developmental  experiments  are  described  in  this 
Appendix.  The  steel  case  of  the  initial  flatpack  design  (used  in  Tests  3  through  5)  was  simply  two 
flat  strips,  without  a  cavity  in  the  base.  The  Kapton/element  sandwich  was  tended  to  the  base, 
and  the  cover  and  base  were  welded  together.  The  lack  of  lateral  confinement  for  the  elements  in 
this  design  led  to  measured  in-plane  strains  as  large  as  4%.  Larger  elements  (~6-mm  square)  were 
also  used  initially,  linuting  the  spatial  resolution  of  the  gage.  Fi^er  modifications  intended  to 
increase  the  frequency  response  of  the  flatpack  were  attempted  with  partial  success  in  Test  6. 
These  are  described  in  Section  4. 


TEST  1 

Figure  D-1  shows  a  sketch  of  the  hardware  we  used  in  otir  first  experimenL  The  4340 
alloy  steel  (Rc  38)  cover  plate  was  1.27  cm  thick.  We  placed  two  sets  of  mangarun  elements  and 
strain-gage  elements  between  the  steel  cover  plate  and  the  steel  cylinder  to  measure  axial  stress  and 
radial  and  circumferential  strain.  One  set  of  gages  was  tended  to  the  cover  plate  and  one  set  was 
tended  to  the  cylinder,  both  located  at  a  2.54-cm  radius.  The  space  between  the  plate  and  cylinder 
was  filled  with  0.33  mm  of  epoxy.  The  explosive  charge  was  a  5.08-cm-diameter,  5.08-cm-tall 
cylinder  of  C4,  at  a  standoff  distance  of  2.54  cm.  The  charge  was  initiated  at  the  center  of  the  top 
surface. 

In  this  experiment,  large  deformations  were  produced  in  both  the  plate  and  the  cylinder, 
and  the  stress  and  strain  gage  elements  were  severely  deformed.  A  S-cm-diameter  depression  was 
made  in  both  the  plate  and  the  cylinder.  The  depression  on  the  top  of  the  plate  was  ateut  0.07-cm 
deep,  and  the  depression  on  the  top  of  the  cylinder  was  ateut  0.03-cm  deep.  The  stress  and  strain 
elements  on  the  plate  were  destroyed  shortly  after  the  arrival  of  the  stress  wave  and  provided  no 
useful  records.  The  gages  on  the  cylinder  survived,  but  the  strain  records  went  off-scale  of  the 
oscilloscopes.  Manganin  gage  records  were  obtained,  but  were  severely  perturbed  by  the  local 
strain. 


TESTS  2  AND  3 

The  configiu-ations  for  Tests  2  and  3  are  shown  in  Figures  D-2  and  D-3,  respec  jvely.  The 
objective  of  the  tests  was  to  obtain  approximate  measurements  of  the  normal  stresses  and  in-plane 
strains  in  these  geometries  to  guide  the  design  of  future  gage  packages.  In  both  tests,  manganin 
and  constantan  elements  were  used.  In  Test  2  the  gage  elements  were  placed  between  the  cylinder 
and  the  cover  plate  (Figure  D-1);  in  Test  3  the  elements  were  packaged  in  flatpacks,  which  were 
cemented  to  the  top  of  the  target  cylinder  (Figure  D-3).  The  target  cylinder  and  cover  plate  were 
made  of  4340  alloy  steel  heat  treated  to  a  harness  of  Rc37  and  ground  flat  and  parallel.  In  both 
tests  the  explosive  charge  was  a  5.08-cm-diameter,  5.08-cm-taIl  cylinder  of  C4  explosive  (165 
grams),  placed  over  the  center  of  the  target  with  a  5.08-cm  standoff. 


D-1 


Figure  D-1 .  Configuration  for  Test  1 . 
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The  manganin  elements  were  Micromeasurements  VM-SS-210AW-048  stress  gages  whose 
active  grid  dimensions  are  approximately  5.3  by  6.3  mm.  The  constantan  elements  were 
Micromeasurements  CEA-()^250UT-350  strain  gages,  whose  active  grid  dimensions  are 
approximately  6.3  by  7.3  mm.  In  Test  2,  sets  of  elements  comprising  one  manganin  stress  gage 
and  two  strain  gages  were  placed  on  the  top  surface  of  the  steel  cylinder  (benea±  the  steel  cover 
plate)  at  2.54-  and  5.08-cm  radii.  In  each  gage  set,  the  gridlines  of  the  stress  gage  and  one  of  the 
strain  gages  were  oriented  in  the  radial  direction;  the  gridlines  of  the  other  strain  gage  were  oriented 
in  the  circumferential  direction.  The  1.27-cm-thick  cover  plate  was  bonded  to  the  cylinder  with  a 
0.2-mm-thick  layer  of  Hysol  epoxy.  The  flatpacks  used  in  Test  3  were  sandwiches  comprising 
two  pieces  of  25-mm-wide  by  1.2-mm-thick  stainless  steel  with  a  0.2-mm-thick  by  12-mm-wide 
epoxy-filled  cavity  between  Aem  in  which  one  stress  gage  and  one  strain  gage  were  mounted  side 
by  side.  The  two  flatpacks  were  bonded  to  the  steel  cylinder  with  epoxy,  with  the  elements 
positioned  at  2.54-  and  5.08-cm  radii  and  oriented  in  ^e  radial  direction. 

The  gage  records  are  plotted  in  Figures  D-4  through  D-7.  The  solid  curves  are  from  Test  2; 
the  dashed  curves  are  from  Test  3.  The  records  from  Test  3  have  been  shifted  by  4  ps  to  maK;h  the 
times  of  arrival  from  Test  2.  The  left  ordinate  scale  is  resistance  change;  the  right  ordinate  scale  is 
approximate  stress  or  strain.  The  latter  are  approximate  because  they  do  not  account  for  the  effects 
of  in-plane  strain  on  the  manganin  elements  or  the  effects  of  normal  stress  on  the  constantan  (strain 
gage)  elements. 

Figure  D-4  shows  that  the  manganin  element  records  from  the  2.54-cm  radius  in  the  two 
tests  are  very  similar  in  shape  and  duration,  but  the  initial  peak  registered  in  the  flatpack  at  the  top 
surface  of  the  target  was  not  registered  under  the  1.27-cm-thick  cover  plate.  (This  difference  is 
qualitatively  consistent  with  the  computational  results  described  below.)  There  is  also  a  small 
amount  of  hysteresis  in  both  elements;  only  slight  additional  resistance  change  occurred  after  the 
time  shown  in  the  plot.  In  general,  these  records  are  encouraging:  the  elements  survived  and 
provided  clean  records  and  there  do  not  appear  to  have  been  large  perturbations  from  in-plane 
strains. 


The  manganin  records  from  the  5.08-cm  radius,  shown  in  Figure  D-5,  suggest  a  much 
different  environment.  The  negative  sense  of  these  records  implies  that  the  effects  of  in-plane 
strain  overwhelmed  the  normal  stress  signals  because  it  is  unreasonable  to  expect  tensile  stresses  at 
these  locations  (compressive  stress  is  positive  in  these  figures).  This  result  suggests  that  elements 
made  of  ytterbium  ,  which  is  much  more  sensitive  to  normal  stress  than  manganin,  could  be 
required  at  this  radius.  Another  observation  is  that  the  two  signals  in  Figure  D-5  are  very  similar  in 
rise  time  and  magnitude,  showing  differences  primarily  in  fi^uency  content  The  record  from 
Test  3  goes  off  scale  due  to  gage  or  lead  failure. 

The  constantan  strain  gage  records  are  shown  in  Figures  D-ti  and  D-7.  Records  from  the 
elements  oriented  in  the  radial  diction  in  Test  2  are  not  shown  because  both  elements  failed  at  the 
time  of  arrival  of  the  first  stress  wave.  All  the  other  records  indicate  that  strains  of  1%  or  greater 
were  registered  at  both  locations  and  in  both  the  radial  and  circumferential  directions.  This  level  of 
strain  would  perturb  the  signals  from  the  manganin  gages  by  the  equivalent  of  a  few  kbar.  The 

strain-gage  records  also  contain  a  strong  oscillatory  component  with  a  period  of  3  to  5  ps.  In  Test 
2  the  oscillations  could  be  associated  with  wave  reverberations  in  the  1.27-cm-thick  cover  plate, 
but  in  Test  3  there  is  no  obvious  source  because  the  cover  plate  was  not  used. 

As  a  check  on  the  consistency  of  the  measurements,  we  can  compare  the  manganin  records 
from  the  5.08-cm  radius  (Figure  D-5),  which  were  apparently  produced  primarily  by  strain,  with 
the  record  from  the  radially  oriented  strain  gage  at  the  same  location  in  Test  3  (Figtire  D-7).  The 

strain  of  approximately  -0.5%  indicated  by  the  strain  gage  at  30  ps  would  produce  a  resistance 
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Figure  D-7.  Strain-gage  records  at  5.08-cm  radius  in  Tests  2  and  3. 

(Tensile  strain  is  positive.) 
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change  in  the  manganin  elements  of  about  -0.7%.  The  manganin  records  indicate  a  resistance 

change  of  about  -0.4%  at  30  (is,  which  leaves  a  difference  of  -0.3%.  This  difference  corresponds 
to  a  stress  of  approximately  1.5  kbar  (0.15  GPa),  which  is  consistent  with  the  calculated 
environment  (see  Appendix  C). 

The  PIEZOR  algorithm  was  used  to  interpret  the  measurements  the  two-element  flatpack 
used  in  Test  3  even  though  the  transverse  strain  was  not  measured.  These  measurements  are 
shown  in  Figure  D-8.  For  this  exercise  we  assumed  that  the  transverse  strain  was  zero.  (In  fact, 
the  transverse  strain  gages  in  Test  5  indicate  that  the  transverse  strain  in  the  flatpack  is  comparable 
to  the  longitudinal  strain,  but  the  transverse  strain  may  be  very  weakly  coupled  to  the  longitudinal 
elements.)  As  shown  in  Figure  D-8,  the  stress  history  computed  with  this  assumption  is  very 
oscillatory.  The  source  of  Aese  oscillations  could  be  rever^ations  in  the  flatpack,although  one 
would  then  expect  them  to  be  more  apparent  in  the  manganin  resistance  history.  The  oscillations 
could  also  be  an  artifact  of  the  assumption  made  in  reducing  the  data  that  the  transverse  strain  is 
zero,  in  which  case  the  effect  of  the  actual  transverse  strain  could  be  to  compensate  for  the 
osciUations  introduced  by  the  longitudinal  strain. 


TEST  4 

The  purpose  of  Test  4  was  to  determine  if  the  large-amplitude  (±2%),  high-frequency 
(200  -  300  kHz)  oscillations  in  the  strain-gage  resistance  histories  in  Tests  2  and  3  could  have  been 
produced  artificially  by  the  signal-conditioning  equipment  whose  linear  response  is  limited  to  about 
100  kHz.  As  shown  in  Figure  D-9,  the  configuration  of  Test  4  was  identical  to  that  of  Test  3 
except  that  only  one  flatpack  gage  was  used.  This  gage  was  the  same  one  used  in  Test  3  at  the 
2.54-cm  radius.  (It  had  survived  Test  3  with  only  surface  damage  and  small  permanent  resistance 
changes  in  the  manganin  and  strain-gage  elements.)  Again  the  gage  was  placed  at  the  2.54-cm 
radius  with  the  elements  oriented  in  the  radial  direction.  The  strain-gage  amplifier,  identical  to 
those  used  with  the  manganin  gages  in  previous  tests,  has  a  linear  fr^uency  limit  of  greater  than  1 
MHz. 


The  unfiltered  strain-gage  record  from  Test  4  shown  in  Figure  D- 10(a)  contains  noise  at  a 
frequency  of  about  3  MHz  in  the  first  10  ps.  This  record  was  filtered  numerically  at  500  kHz  to 
obtain  a  more  direa  comparison  with  the  strain  record  from  Test  3.  As  shown  in  Figures  D-10(b) 
and  (c),  the  filtered  strain  gage  record  from  Test  4  has  a  different  character  than  the  record  obtained 
in  Test  3.  Unfonunately,  not  only  is  the  frequency  content  different  but  the  overall  amplitude  is 
also  larger.  Although  a  nonlinear  filter  applied  to  the  Test  4  record  might  produce  a  closer  match 
to  the  Test  3  record,  we  believe  the  element  in  Test  4  experienced  a  different  strain  environment,  so 
the  different  response  cannot  be  attributed  solely  to  different  signal  conditioning.  The  mechanical 
properties  of  the  flatpack  may  have  been  affected  by  its  use  in  an  earlier  test.  However,  the  Test  4 
record  clearly  has  significant  frequency  content  at  about  150  kHz  and  greater.  Thus,  the  high- 
frequency  signal  conditioning  for  strain  gages  was  used  in  subsequent  experiments. 


RESIDUAL  DEFORMATION  PROFILES 

The  purposes  of  measuring  the  deformed  shapes  of  the  steel  targets  are  to  help  deteimine 
the  strain  environment  that  the  gages  experienced  in  ^e  tests  and  to  help  determine  the  accuracy  of 
the  computational  predictions  presented  in  Appendix  C.  The  deformation  profiles  were  measured 
by  traversing  the  surface  of  the  test  article  along  two  orthogonal  diameters  with  a  dial  indicator 
mounted  in  ^e  head  of  a  milling  machine.  The  four  elevation  measurements  taken  at  each  radius 
were  then  averaged.  For  Test  2  we  measured  the  shapes  of  the  top  surface  of  the  plate,  the  bottom 
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Figure  D-8.  PIEZOR  interpretation  of  Test  3  flatpack  records  for  zero 
transverse  strain. 
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Figure  0-9.  Configuration  for  Test  4. 


surface  of  the  plate,  and  the  top  surface  of  the  cylinder.  For  Tests  3  and  4,  we  measured  the  shape 
of  the  top  surface  of  the  cylinder  after  Test  4.  Tlius,  this  shape  contains  the  accumulated 
permanent  deformation  from  Tests  3  and  4.  In  the  plots  shown  below,  these  accumulated 
deformations  were  halved  to  approximate  those  from  a  single  loading. 

Figure  D-1 1(a)  shows  the  deformed  profile  of  the  top  surface  of  the  plate  after  Test  2.  The 
depth  of  the  depression  at  the  center  of  the  plate  is  about  0.016  cm  (6.3  mil).  The  computed 

profiles  are  from  the  last  time  step  in  the  computations  (50  jis).  Both  computations  (for  a  no-slip 
interface  and  a  free-slip  interface)  predict  a  central  depression  of  about  the  same  diameter  as 
measured,  but  the  predicted  depth  of  the  depression  is  much  larger  than  measured.  The 
computations  also  predict  much  more  deformation  outside  the  central  depression  than  is  measured. 

Figure  D-1 1(b)  shows  the  deformed  shapes  of  the  bottom  of  the  plate  and  the  top  of  the 
cylinder  after  Test  2.  By  spacing  these  profiles  0.020  cm  (8  mil)  apart  at  the  outer  radius,  we  have 
an  approximate  view  of  the  residual  shape  of  the  epoxy  layer.  The  epoxy  layer  apparently 
compressed  quite  a  bit  near  the  center,  which  must  have  significantly  affected  the  deformatiem  of 
the  plate  and  cylinder.  Also,  the  gradient  in  residual  thickness  of  the  epoxy  layer  suggests  that  the 
plate  deformation  produced  radial  flow  in  the  epoxy  layer,  which  could  have  ^ected  the  strain 
gages  and  manganin  elements  on  the  surface  of  the  cylinder.  Both  computational  simulations 
overpiedict  the  deformation  of  the  cylinder,  the  computations  make  no  prediction  of  residual  epoxy 
layer  thickness  because  they  both  modeled  the  epoxy  layer  as  an  interface. 

Figure  D- 1 1(c)  shows  the  deformed  shape  of  the  top  of  the  cylinder  after  Tests  3  and  4, 
with  the  measured  deformations  halved  to  approximate  the  result  of  a  single  loading.  Again,  the 
computational  simulation  (with  the  no-slip  interface)  overpredicts  the  depth  of  the  central 
depression  and  predicts  a  defamed  shape  outside  the  central  depression  that  is  unlike  the  observed 
profile. 


The  deformed  profile  measurements  from  Tests  2  to  4  indicate  that  significant  residual 
deformation  occurs  with  a  4-cm  radius,  so  it  would  be  reasonable  to  expect  large  strains  in  that 
portion  of  the  target.  (In  future  experiments,  we  will  scribe  the  targets  with  equally  spaced  lines 
before  the  test  so  that  we  can  make  posttest  measurements  of  residual  strain  in  addition  to  residual 
surface  shape.)  However,  the  computational  simulations  oveipredict  the  deformation,  suggesting 
that  the  computed  strains  are  also  greater  than  those  experienced  by  the  target.  Thus,  the  strain- 
gage  signals,  which  greatly  exceed  the  predicted  strains,  apparently  do  not  represent  the  actual 
strains  in  the  target.  RecaU  that  an  elastic  perfectly  plastic  constitutive  model  was  used  for  the 
cylinder  and  cover  plate.  Not  allowing  h^ening  could  lead  to  an  overestimation  of  the  residual 
strains. 


ASSESSMENT  OF  THE  STRAIN-GAGE  RECORDS 

The  strains  predicted  by  L2D  computations  (Appendix  C)  are  about  an  order  of  magnitude 
smaller  than  those  indicated  by  the  strain-gage  records  from  Tests  2  and  3.  However,  as  indicated 
in  the  previous  subsection,  it  does  not  seem  plausible  that  the  target  strains  were  10  times  larger 
than  predicted  because  the  measured  residual  deformations  were  considerably  smaller  than  the  L2D 
predictions.  Thus,  the  strain-gage  signals  were  not  induced  directly  by  target  strains. 

The  magnitudes  of  normal  stress  oscillations  that  would  be  required  to  produce  the 
observed  strain-gage  oscillations  are  very  large.  In  the  case  of  the  circumferential  strain  gage  at  the 
2.54-cm  radius  in  Test  2,  the  stress  oscillations  would  have  to  have  been  25  kbar  (2.5  GPa). 

Thus,  the  strain-gage  signals  were  not  induced  by  stress  oscillations. 
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Figure  D-l  1 .  Residual  deformation  profiles  from  Tests  2, 3,  and  4. 
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Oscillations  of  the  same  frequency  as  those  in  the  strain-gage  signals  also  occurred  in  the 
manganin-gage  signals  in  Test  3.  Furthermore,  the  amplitudes  of  the  oscillations  in  the  manganin- 
gage  signals  are  consistent  with  the  assumption  that  the  strain-gage  signals  represent  elastic  strain 
in  the  manganin  element  Thus,  it  appears  that  the  strain  gages  are  recording  the  actual  clement 
environment 

The  frequency  of  the  oscillations  in  the  strain-gage  signals  from  Test  2  matches  the 
frequency  of  stress  oscillations  in  the  epoxy  layer  pre^cied  by  a  one-dimensional  stress-wave 
computation  (not  shown).  Also,  the  oscillations  in  the  strain-gage  signals  from  the  flatpacks  in 
Test  3  approximately  match  the  calculated  flexural  frequency  of  the  flatpack  cover  supported  by  an 
elastic  layer  of  epoxy.  Thus,  it  appears  that  package  response  (the  plate  on  the  ejwxy  layer  or  the 
flatpack  on  the  surface)  is  inducing  strains  in  the  elements  that  are  much  different  in  character  and 
magnitude  from  the  local  target  strains. 

In  summary,  the  evidence  indicates  that  the  strain  gages  in  Tests  2  and  3  provide  records  of 
element  strains  induced  by  the  response  of  the  element  packages.  That  is,  the  strain-gage  reewds 
represent  strains  in  the  epoxy  layer  in  Test  2  and  strains  inside  the  flatpacks  in  Test  3.  That  these 
strain  histories  are  so  different  (much  larger  amplitude)  fhun  the  computational  predicdcxis  of  target 
strain  implies  that  the  elements  need  to  be  packaged  or  positioned  differently  to  inhibit  these 
strains. 


TEST  5 

The  objectives  of  Test  5  were  to 

( 1 )  Determine  the  relative  magnitudes  of  strain  in  the  longitudinal  and  transverse 
directions  in  flatpack  gages. 

(2)  Find  the  dependence  of  these  strains  on  the  nature  of  the  attachment  of  the 
flatpack  to  the  target  and  the  orientation  of  the  flatpack  to  the  explosive  charge. 

The  25-cm-diameter  steel  target  and  the  5-cm-diameter  C4  explosive  charge  with  a  5-cm 
standoff  were  the  same  as  in  the  previous  tests.  Four  flatpack  gages  were  used,  e^h  with  a  pair  of 
orthogonal  strain-gage  elements.  Figure  D- 12  shows  the  flatpack  gage  configuration.  Two 
flatpacks  were  weld^  to  the  target;  two  were  bonded  with  epoxy.  Two  were  oriented  radially  to 
the  charge;  two  were  oriented  tangentially.  The  center  of  the  sensing  area  of  each  flatpack  was  at  a 
2.75-cm  r^us  from  the  axis  of  the  explosive  charge.  (The  slightly  larger  than  2.54-^  spacing 
was  required  to  allow  room  for  welding.) 

The  strain-gage  records  are  plotted  in  Hgures  D-13  through  D-16.  In  these  plots,  the 
recorded  resistance  histories  have  b^n  digitally  filtered  at  1  MHz  and  converted  to  strain  using  the 
manufacturer’s  gage  calibration  without  compensating  for  normal  stress.  In  general,  the  strain  in 
the  transverse  dvection  is  comparable  in  magnitude  to  that  in  the  longitudinal  direction.  Also,  the 
welded  flatpacks  had  slightly  lower  strains  ^an  those  bonded  with  epoxy.  Radially  oriented  and 
tangentially  oriented  flatpacks  had  about  the  same  strain  levels. 

From  these  results  we  conclude  that 

(1)  Both  longitudinal  and  transverse  strains  are  significant.  Thus,  both  strains 
need  to  measured  to  account  for  the  effects  of  in-plane  strain  on  manganin 
and  ytterbium  stress-sensing  elements. 


D-16 


(a)  Top  view  of  gage  arrangement 


(b)  Flatpack  design 
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Figure  D-1 2.  Flatpack  gage  configuration  for  Test  5. 
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Figure  D-13.  Test  5  strains  from  the  welded  flatpack  in  the  radial  orientation. 
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(2)  The  strains  in  welded  flatpacks  are  lower  than  those  in  epoxied  flatpacks,  but 
they  are  still  large  enough  to  significantly  affect  stress-sensing  elements. 

(3)  The  orientation  of  the  flatpack  relative  to  the  explosive  charge  is  not 
important 
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APPENDIX  E 


SIGNAL  CONDITIONING 


The  pulsed  bridge,  shown  schematically  in  Figure  E-1,  is  a  high-voltage  version  of  an 
ordinary  Wheatstone  resistance  bridge  in  which  a  change  in  resistance  is  determined  by  measuring 
the  voltage  due  to  bridge  imbalance.  The  higher  the  voltage  applied  to  the  bridge,  the  higher  the 
imbalance  signal.  The  design  of  this  circuit  is  described  next;  the  relation  between  output  voltage 
and  gage  resistance  is  discussed  later. 

To  avoid  destroying  the  piezoresistance  gage  element  by  the  high  current  and  resultant 
heating,  the  duration  of  bridge  voltage  application  is  limited  to  Ae  time  needed  for  the 
measurement.  This  can  be  done  by  discharging  a  capacitor  (Ci)  with  a  triggered  silicon-controlled 
rectifier  (SCR)  switch  to  apply  the  voltage  and  a  crowbar  circuit  to  turn  off  the  voltage,  as  done  in 
SRI  power  supplies.  Alternatively,  one  can  use  a  field-effect  transistor  (FET)  to  control  voltage, 
as  done  in  SAIC*  power  supplies. 

The  apparent  simplicity  of  the  pulsed  bridge  circuit  is  deceptive.  This  is  an  AC  circuit, 
designed  to  measure  over  a  wide  range  of  frequencies  the  change  in  resistance  of  a  gage  element 
located  at  the  end  of  a  long  or  short  cable.  Funhermore,  the  gage  element  is  in  a  generally  noisy 
environment,  and  the  system  must  be  optimized  ior  noise  rejection. 

Referring  to  Rgure  E-1,  Ri  and  R2  are  10  ohm  resistors  used  to  limit  current  through  the 
crowbar  SCR.  They  should  be  at  least  5  watt  resistors;  10  watt  would  be  safer.  They  arc  not  used 
or  needed  in  the  SAIC  pulsed  power  supply,  which  uses  a  power  FET  rather  than  an  SCR  crowbar 
to  turn  the  power  off. 

For  signals  traveling  fiom  the  gage  to  the  power  supply,  the  capacitor  Ci  appears  as  a  short 
circuit.  To  properly  terminate  the  gage  cable,  the  sum  of  Ri  +  R2  +  Rs  should  be  equal  to  the 
characteristic  impedance  of  the  gage  cable  (within  10%).  IdeaUy,  the  gage  cable  should  be  a  good 
quality  balanced  twin-ax  such  as  RG22B.  RG22B  cable  has  a  characteristic  impedance  of  95 
ohms.  With  Ri  =  R2  =  10  ohms,  one  would  choose  a  value  of  75  ohms  for  Rg.  With  the  SAIC 
supply,  Ri  +  R2  =  0,  and  Rg  would  be  95  ohms. 

Gage  element  resistance  is  a  design  option.  Since  it  is  the  relative  resistance  change  that  is 
related  to  the  magnitude  of  the  induced  stress,  the  magnitude  of  the  resistance  change  at  a  given 
stress  will  be  simply  proportional  to  the  initial  resistance  of  the  gage  element.  If  the  current- 
hwdling  capacity  were  independent  of  gage  resistance,  one  would  expect  larger  signals  fiom 
higher  initial  resistance  gage  elements.  However,  if  the  area  of  the  gage  element  is  held  constant, 
the  current-handling  capacity  of  the  gage  decreases  proportionately  as  the  resistance  increases. 

Gage  cable  resistance  should  be  as  low  as  practical;  the  sensitivity  of  the  measurement 
system  is  inversely  related  to  the  gage  cable  resistance.  C^ble  quality  is  important;  it  is  imperative 
for  accurate  data  i^uction  to  c^ture  the  peak  resistance  changes  by  using  high  quality  cable  and/or 
shorter  runs  to  the  instrumentation. 


’Sources  for  pulsed-power  supplies  are  listed  at  the  end  of  this  section. 
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The  sum  of  R3  +  R5  +  R4  is  generally  in  the  range  of  700  to  1500  ohms  for  most 
measurements.  This  range  of  values  is  a  compromise:  iugh  enough  to  permit  use  of  inexpensive 
one-watt  ten-turn  potentiometers  but  not  so  high  that  the  high-frequency  response  of  the  system  is 
degraded  by  stray  capacitance.  Two  hundred  ohms  is  a  go^  value  for  Rs,  the  ten-turn  pot.  It  is 
necessary  to  choose  values  of  R3  and  R4  compatible  with  the  ratio  of  Rg  to  the  sum  of  gage  and 
cable  resistances.  If  the  sum  of  gage  and  cable  resistance  is  about  20  ohms  and  Rg  is  75  ohms  (for 
impedance-matching  reasons),  then  R3  and  R4  could  be  750  and  200  ohms,  respectively.  If  the 
sum  of  gage  and  cable  resistance  is  1()0  ohms  and  I^  is  75  ohms,  R3  and  R4  could  be  500  and  600 
ohms,  respectively.  We  generally  use  two-watt  resistors  for  R3  and  R4.  As  mentioned,  R5  is  a 
one- watt  ten-turn  pot. 

R6  is  zero;  we  don't  use  a  resistor  here  at  all. 

R7  is  the  common  mode  balance  resistor: 

I/R7  =  l/(Ri  +  R2  +  R3  +  Rx)  +  1/[R4  +  (Rs  -  Rx)] 

where  Rx  (the  part  of  the  pot  resistance  on  the  R3  side  of  the  slider)  can  be  calculated  from  the 
relationstdp: 


(R3  +  Rx)/(R3  +  Rs  +  R4)  =  R8/(R8  +  R  gage  +  R  cable). 

The  purpose  of  R7  is  to  balance  the  common  mode  signals  on  the  gage  cable  with  respect  to 
the  inputs  to  the  differential  amplifier.  Noise  generated  in  the  test  configuration  will  be  coupled 
from  the  shield  to  the  gage  cable  as  common  mode  noise.  If  R7  is  zero,  one  side  of  the  cable  will 
be  coupled  direcdy  to  the  input  impedance  of  one  side  of  the  amplifier.  The  other  side  of  the  cable 
will  be  coupled  to  the  other  side  of  the  amplifier  through  the  parcel  combination  of  (Ri  +  R2  + 

R3  +  Rx)  and  0^4  +  Rs  -  Rx)-  If  the  input  impedance  of  the  amplifier  is  very  high  (megohms), 

R7  may  be  omitted.  If  the  input  impedance  is  reduced  to  the  Kilohm  range,  e.g.,  through  use  of 
anenuators,  then  the  common  mode  balance  obtained  by  choosing  an  appropriate  value  for  R7  can 
be  critical  in  a  noisy  environment.  Note  that  the  value  of  R7  was  about  300  ohms  and  served  to 
correct  a  3%  imbalance.  A  10%  error  in  the  value  of  R7  was  acceptable,  corresponding  to  an 
imbalance  of  0.3%.  R7  may  be  a  1/8-watt  resistor. 

The  purpose  of  the  attenuator  is  to  protect  the  amplifier  from  both  common-mode  and 
normal-mode  overvoltages.  Because  of  cable  and  system  capacitance  to  ground,  a  decaying 
common  nnode  voltage  appears  at  the  amplifier  inputs  when  the  bridge  circuit  is  pulsed.  The  initial 
level  of  this  common  mcde  voltage  is  determined  by  the  ratio  of  the  gage  and  cable  resistance  to  Rg 
and  can  exceed  the  common  mode  voltage  handling  capabilities  of  the  differential  amplifier. 
Additiotutlly,  a  large  normal  mode  spike  can  be  genoated  when  the  supply  is  pulsed  because  of  the 
AC  imbalance  of  the  bridge.  The  attenuator  reduces  normal-  and  common-mc^e  voltages  to  levels 
that  can  be  handled  by  the  amplifier.  We  typically  use  a  10:1  attenuator.  Note  that  the  attenuators 
are  made  up  in  matched  pairs.  Common  mode  rejection  can  be  seriously  compromised  by  mixing 
up  the  originally  matched  pairs.  Our  attenuator  recipe  is: 

R9  =  Ri2  These  resistors  are  nominally  9  k  ohms  and  are  matched  to  within 
0.05%. 

RlO  =  Rl  1  These  resistors  are  nominally  1  k  ohm  and  are  matched  to  within 
0.05%. 

1/8-watt  resistors  are  adequate  for  the  attenuators. 
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The  attenuators  are  carefully  matched  so  that  the  common  mode  noise  will  be  attenuated 
equally  at  both  amplifier  inputs.  In  principle,  matched  capacitors  should  be  added  to  the 
attenuators  to  improve  their  high  frequency  characteristics.  However,  after  conducting  tests  at  SRI 
with  both  capacitor-compensated  and  uncompensated  attenuators,  we  concluded  that  the  capacitors 
made  no  detectable  difference,  given  the  configuration  with  short  cables  between  the  power 
supplies  and  amplifiers.  Note  that  R7  must  be  included  in  an  accurate  calculation  of  the  magnitude 
of  the  attenuation.  With  R7  about  300  ohms,  the  attenuation  will  be  increased  by  about  1.5%. 

The  source  voltage  Vs  applied  to  the  bridge  circuit  by  the  capacitor  Ci  is  given  by 


Vs  =  VsoC-‘A 


where  Vso  is  the  initial  capacitor  voltage,  and  t  is  the  time  constant  given  by 


T  =  Ri  +  R2  +  Rab 


where  Rab  is  the  resistance  between  A  and  B  of  the  bridge,  given  by 

[Ratten  (R:^+R4+R5)^^(R3+Rx)(R4-^R5-Rx)]  (R«+Rg-<-Rc)+Rr  (R^.■^Rc)(R3^^R4+R.s) 
~  Ratten  (R8+R3+R4+R5+Rg+Rc)+(R8+R3+Rx)(Rg+Rc+R4+R5-Rx) 

where  Rattan  =  R6  +  R9  +  R10  +  R11 +  R12  +  R7 
Rg  =  Gage  resistance 

Rc  =  Cable  resistance 

The  time  constant  t  is  determined  at  the  arrival  time  (t  =  tg)  of  the  shock  wave  by  assuming  only 
the  initial  value  of  Rg-  We  have  found  that  changes  in  Rg  arising  from  the  shock  wave  have  no 

appreciable  effect  on  Vj  through  changes  in  x,  at  times  greater  than  to- 

The  bridge  voltage  is  Vcd  is  related  to  the  differential  amplifier  output  voltage  Ve  by 

,,  .  Attenuation  factor 

VcD  =  VEX^-jTjj^^ 


The  attenuation  factor  is  given  by 


Attenuation  factor  — — iP— 

fatten 

but  attenuation  factor  and  gain  are  usually  determined  by  measurement. 

We  are  interested  in  the  effect  of  changes  in  Rg  arising  from  the  shock  wave.  Because  we 
measure  Ve  and  determine  Vcd,  we  wish  to  determine  Rg  from  Vcd,  so  we  can  then  find  the 
stress  on  the  piezoresistance  gage.  The  expression  for  Rg  is  given  by 
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Rg  =  -Rc  +  {  Vs  Ratten  (R4  +  Rs  “  Rx)  Rs  -  VcD  [-2(R3+Rx)  Ratten  (R4+R5-RX) 

+  (R4+R5-Rx+Ratten)  (Rs+Rx+Rs+Ratlen)  *  (R3+R4+R5+R1+R2) 

-  (R3+Rx)2  (R4+R5-Rx+Ratten)  -  (Rg+Rs+Rx+Ratlen)  (R^+Rs-Rx)^ 

~  Ratten^(R3+R4+R5+Rl+R2)]  }"*" 


{  Vs  Ratten  (Rs+Rx)  +  VcD  [(Rg+Rs+Rx+Ratten)  (R3+R4+R5+R1+R2) 

-  (R3+Rx)2]  } 

In  general,  the  bridge  will  not  be  perfecdy  balanced,  because  of  limited  resolution  of  the 
potentiometer,  because  of  only  an  estimate  for  the  expected  change  in  Rg,  and  because  of  ambient 
heating  resistance  changes  in  cable  and  gage.  The  offset  or  reference  voltage  is  referred  to  as 
VcDO.  observed  at  t  =  to.  We  can  now  obtain  the  change  in  gage  resistance  ARq  with  respect  to  its 
reference  value  by  the  following  expression: 


aRg  =  Rg  (VcD(t).  Vs(t),  Ri)  -  Rg  (Vcdo.  Vs(to).  RO 

assuming  no  further  resistance  heating  changes  occur  in  cable  or  gage  (or  elsewhere  in  the  circuit). 
The  relative  change  in  gage  resistance  is  then  given  by 

ARg 

Rgo 

where  Rgo  is  the  gage  resistance  at  t  =  to.  If  gage  heating  effects  are  negligible,  Rgo  will  be  the 
ambient  temperature  value  of  Rg- 


SOURCES  FOR  PULSE  POWER  SUPPLIES 

1 .  SAIC  -  Science  Applications  International  Corporation,  3351  South  Highland  Drive,  Suite 
206,  Las  Vegas,  NV  89109,  Model  SG-133A. 

2.  Dynasen,  Inc.,  20  Arnold  Place,  Goleta,  CA  931 17  (805)  964-4410,  Model  CKl-50-300. 
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APPENDIX  F 


PRETEST  CALCULATIONS  OF  THE  PLATE-IMPACT  EXPERIMENT 


In  the  plate-impact  experiment  described  in  Section  5,  a  25.4-mm-thick  PMMA 
flyer  plate  inpacted  a  2S.4-mm-thick  buffer  plate.  Flatpack  gages  were  placed  at  the 
interface  between  the  buffer  plate  and  a  steel  target  plate.  In  the  design  of  this  experiment, 
the  concern  was  that  the  flatpack  gages  might  ui^er-iegister  the  interface  stress  b^use 
they  are  more  compliant  than  the  solid  steel  target.  In  this  appendix,  we  present  the  results 
of  one-,  two-,  and  three-dimensional  pretest  calculations  performed  for  the  specific  purpose 
of  predicting  the  response  of  the  flatpacks. 

In  these  calculations,  the  steel  was  modeled  as  an  elastic-plastic  material.  The 
piezoresistance  elements,  the  Kapton  insulation,  and  the  epoxy  inside  the  flatpacks  were 
modeled  as  a  unifoim  layer  with  the  elastic-plastic  properties  of  the  epoxy.  The  PMMA 
was  linearly  elastic.  Table  F-1  gives  the  properties  of  these  materials.  All  the  interfaces 
between  materials  were  bonded.  The  impact  velocity  in  the  calculations  was  SI  m/s  (higher 
than  actually  achieved  in  the  experiment).  The  one-dimensional  calculations  were 
performed  with  SRI  PUFF.*  The  two-  and  three-dimensional  calculations  were  performed 
with  DYNA3D.** 


Table  F-1.  Material  Properties  for  Pretest  Calculations 


Material 

Density 

(g/cm3) 

Young's  Modulus 
(GPa) 

Poisson’s 

Ratio 

Yield  Strength 
(GPa) 

Steel 

7.85 

200 

0.3 

1.24 

1.0 

PMMA 

1.184 

6.73* 

0.345 

1. 

1.50 

3.45"* 

0.35 

0.069"* 

0.1 

*The  elastic  behavior  and  the  value  of  Young's  modulus  represent  the  estimated  behavior  for  the 
strain  rates  involved.  The  static  compressive  strength  is  about  0.083  GPa. 

**The  manufacturer's  recommended  values  are  3.17  and  0.0B3  GPa. 


*L.  Seaman,  "SRI  PUFFS  Computer  Program  for  One-Dimensional  Stress  Wave  Propagation,"  SRI 
International  Final  Report  fen  U.S.  Army  Ballistic  Research  Laboratmy,  Contract  DAAKl  1-77-C-0083, 
Vol.  II  (August  1978). 

**J.  Hallquist  and  D.  Benson.  "DYNA3D  User's  Manual  (Nonlinear  Dynamic  Analysis  of  Structures  in 
Three  Dimensions),"  Lawrence  Livomore  National  Laboratory.  Report  UCID-19S92,  Rev.  2  (March 
1986). 
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Figure  F-1  shows  the  results  of  a  pair  of  one-dimensional  calculations.  The  first 
calculation  did  not  include  the  flatpack  gage.  Figure  F-1  (a)  shows  the  computed  interface 
stress  for  this  case.  In  the  second  calculation  a  one-dimensional  flatpack  was  modeled  with 
a  0.94-mm-thick  steel  layer  and  a  0.36-mm-thick  epoxy  layer  between  the  PMMA  buffer 
plate  and  the  steel  target.  Figure  F-l(b)  shows  that  an  overshoot  and  dampal  ringing  of  the 
stress  inside  the  gage  are  predicted,  but  that  the  nominal  equilibrium  stress  inside  the  gage 
is  essentially  equal  to  the  interface  stress. 

Figure  F-2  shows  the  results  of  a  two-dimensional  calculation  in  which  the  flatpack 
steel  case  and  epmxy-filled  cavity  were  included.  The  plane  of  this  calculation  included  the 
normal  and  transverse  axes  of  the  flatpack,  that  is,  it  is  a  plane  strain  calculation  with  the 
longitudinal  strain  in  the  flatpack  prolubited.  The  three  cells  for  which  the  output  is  plotted 
are  in  the  epoxy-filled  cavity  of  the  flatpack  gage.  One  cell  (A)  borders  the  center  of  the 
sensing  area,  one  cell  (B)  borders  the  outer  ^gc  of  the  sensing  area,  and  one  point  (C) 
borders  the  ^ge  of  the  cavity.  The  stress  plots  show  that  within  the  sensing  area  there  is 
an  overshoot  followed  by  oscilladons  about  the  equilibrium  stress  level,  but  that  the 
equilibrium  stress  level  is  essentially  the  same  as  Ae  one-dimensional  interface  stress. 

Near  the  edge  of  the  cavity  the  stress  inside  the  gage  is  much  lower.  The  strain  plots  show 
that  the  transverse  strain  is  less  that  0.002  within  the  sensing  area  and  is  only  slightly  larger 
near  the  edge  of  the  cavity. 

The  results  of  a  three-dimensional  calculation  are  shown  in  Figiue  F-3.  Both  cells 
for  which  the  output  is  plotted  are  in  the  epoxy-filled  cavity  of  the  flatpack.  One  cell  (A)  is 
at  the  center  of  the  sensing  area;  the  other  cell  (B)  is  27.5  mm  from  the  sensing  area 
(toward  the  outer  edge  of  the  target).  Both  stresses  show  an  overshoot  and  oscillations. 
The  normal  stress  at  A  is  slightly  higher  than  the  equilibrium  stress  predicted  by  the  one- 
and  two-dimensional  calculations;  the  normal  stress  at  B  is  slightly  lower  than  those 
predictions.  The  longitudinal  and  transverse  strains  are  less  than  0.0005  except  for  one 
excursion  of  longinuhnal  strain  at  B  to  about  0.001. 

In  summary,  the  pretest  calculations  do  not  suggest  that  significant  under¬ 
registration  of  stress  by  the  flatpacks  should  occur  in  the  plate-impact  experiment 
Additional  analysis  is  necessary  to  fully  understand  the  results  of  the  experiment. 
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(b)  Normal  stress  inside  the  flatpack 

RA-7422-73 

Figure  F-1 .  One-dimensional  finite-element  pretest  calculations 
of  the  plate-impact  experiment. 


STRESS  (kbar) 


A  «  Flatpack  Centerline 
B  *  Element  Outer  Edge 
C  *  Cavity  Edge 


10  2 
TIME  (us) 

(a)  Normal  stress  inside  the  flatpack 


Figure  F-2.  Two-dimensional  finite-element  pretest  calculation 
of  the  plate-impact  experiment. 
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(b)  Transverse  strain  inside  the  flatpack 
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Rgure  F-2.  Two-dimensional  finite-element  pretest  calculation 
of  the  plate-impact  experiment.  (Concluded) 
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Figure  F*3.  Three-dimensional  finite-element  pretest  calculation 
of  the  plate-impact  experiment. 
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Figure  F-3.  Three-dimensional  finite-element  pretest  calculation 
of  the  plate-impact  experiment.  (Continued) 
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(c)  Transverse  strain  inside  the  flatpack 


Figure  F-3.  Three-dimensional  finite-element  pretest  calculation 
of  the  plate-impact  experiment.  (Concluded) 
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